! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! The Reaction Rates File
! 
! Generated by KPP-2.2.3 symbolic chemistry Kinetics PreProcessor
!       (http://www.cs.vt.edu/~asandu/Software/KPP)
! KPP is distributed under GPL, the general public licence
!       (http://www.gnu.org/copyleft/gpl.html)
! (C) 1995-1997, V. Damian & A. Sandu, CGRER, Univ. Iowa
! (C) 1997-2005, A. Sandu, Michigan Tech, Virginia Tech
!     With important contributions from:
!        M. Damian, Villanova University, USA
!        R. Sander, Max-Planck Institute for Chemistry, Mainz, Germany
! 
! File                 : aromatics_kpp_Rates.f90
! Time                 : Fri Jun 12 17:21:21 2020
! Working directory    : /n/home08/kbates/Aromatics/MECCA
! Equation file        : aromatics_kpp.kpp
! Output root filename : aromatics_kpp
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



MODULE aromatics_kpp_Rates

  USE aromatics_kpp_Parameters
  USE aromatics_kpp_Global
  IMPLICIT NONE

CONTAINS



! Begin Rate Law Functions from KPP_HOME/util/UserRateLaws

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  User-defined Rate Law functions
!  Note: the default argument type for rate laws, as read from the equations file, is single precision
!        but all the internal calculations are performed in double precision
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

!~~~>  Arrhenius
   REAL(kind=dp) FUNCTION ARR( A0,B0,C0 )
      REAL A0,B0,C0      
      ARR =  DBLE(A0) * EXP(-DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
   END FUNCTION ARR        

   REAL(kind=dp) FUNCTION TUN( A0,B0,C0 )
      REAL A0,B0,C0      
      TUN =  DBLE(A0) * EXP(DBLE(B0)/TEMP) * EXP(DBLE(C0)/TEMP**3)
   END FUNCTION TUN

   REAL(kind=dp) FUNCTION ISO1( A0,B0,C0,D0,E0,F0,G0 )
      REAL A0,B0,C0,D0,E0,F0,G0
      REAL(kind=dp) K0, K1, K2
      K0 = DBLE(D0)*EXP(DBLE(E0)/TEMP)*EXP(1.E8/TEMP**3)
      K1 = DBLE(F0)*EXP(DBLE(G0)/TEMP)
      K2 = DBLE(C0)*K0/(K0+K1)
      ISO1 =  DBLE(A0) * EXP(DBLE(B0)/TEMP) * (1.-K2)
   END FUNCTION ISO1

   REAL(kind=dp) FUNCTION ISO2( A0,B0,C0,D0,E0,F0,G0 )
      REAL A0,B0,C0,D0,E0,F0,G0
      REAL(kind=dp) K0, K1, K2
      K0 = DBLE(D0)*EXP(DBLE(E0)/TEMP)*EXP(1.E8/TEMP**3)
      K1 = DBLE(F0)*EXP(DBLE(G0)/TEMP)
      K2 = DBLE(C0)*K0/(K0+K1)
      ISO2 =  DBLE(A0) * EXP(DBLE(B0)/TEMP) * K2
   END FUNCTION ISO2

   REAL(kind=dp) FUNCTION ARRPLUS( A0,B0,C0,D0,E0 )
      REAL A0,B0,C0,D0,E0 
      REAL(kind=dp) K0     
      K0 = DBLE(A0) * (DBLE(D0)+(TEMP*DBLE(E0)))
      ARRPLUS =  K0 * EXP(-DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
      IF (ARRPLUS < 0.0) ARRPLUS = 0.0
   END FUNCTION ARRPLUS

   REAL(kind=dp) FUNCTION TUNPLUS( A0,B0,C0,D0,E0 )
      REAL A0,B0,C0,D0,E0
      REAL(kind=dp) K0
      K0 = DBLE(A0) * (DBLE(D0)+(TEMP*DBLE(E0)))
      TUNPLUS = K0 * EXP(DBLE(B0)/TEMP) * EXP(DBLE(C0)/TEMP**3)
      IF (TUNPLUS < 0.0) TUNPLUS = 0.0
   END FUNCTION TUNPLUS

!~~~> Simplified Arrhenius, with two arguments
!~~~> Note: The argument B0 has a changed sign when compared to ARR
   REAL(kind=dp) FUNCTION ARR2( A0,B0 )
      REAL A0,B0           
      ARR2 =  DBLE(A0) * EXP( DBLE(B0)/TEMP )              
   END FUNCTION ARR2          

   REAL(kind=dp) FUNCTION EP2(A0,C0,A2,C2,A3,C3)
      REAL A0,C0,A2,C2,A3,C3
      REAL(kind=dp) K0,K2,K3            
      K0 = DBLE(A0) * EXP(-DBLE(C0)/TEMP)
      K2 = DBLE(A2) * EXP(-DBLE(C2)/TEMP)
      K3 = DBLE(A3) * EXP(-DBLE(C3)/TEMP)
      K3 = K3*CFACTOR*1.0E6_dp
      EP2 = K0 + K3/(1.0_dp+K3/K2 )
   END FUNCTION EP2

   REAL(kind=dp) FUNCTION EP3(A1,C1,A2,C2) 
      REAL A1, C1, A2, C2
      REAL(kind=dp) K1, K2      
      K1 = DBLE(A1) * EXP(-DBLE(C1)/TEMP)
      K2 = DBLE(A2) * EXP(-DBLE(C2)/TEMP)
      EP3 = K1 + K2*(1.0E6_dp*CFACTOR)
   END FUNCTION EP3 

   REAL(kind=dp) FUNCTION EPO(A1,E1,M1)
      REAL A1, E1, M1
      REAL(kind=dp) K1      
      K1 = 1.0_dp/(DBLE(M1) * CFACTOR + 1.0_dp)
      EPO = DBLE(A1) * EXP(DBLE(E1)/TEMP) *  K1
   END FUNCTION EPO

   REAL(kind=dp) FUNCTION KCO(A1,M1)
      REAL A1, M1
      KCO = DBLE(A1) * (1.0_dp + (CFACTOR / DBLE(M1)))
   END FUNCTION KCO

   REAL(kind=dp) FUNCTION FALL ( A0,B0,C0,A1,B1,C1,CF)
      REAL A0,B0,C0,A1,B1,C1,CF
      REAL(kind=dp) K0, K1     
      K0 = DBLE(A0) * EXP(-DBLE(B0)/TEMP)* (TEMP/300.0_dp)**DBLE(C0)
      K1 = DBLE(A1) * EXP(-DBLE(B1)/TEMP)* (TEMP/300.0_dp)**DBLE(C1)
      K0 = K0*CFACTOR*1.0E6_dp
      K1 = K0/K1
      FALL = (K0/(1.0_dp+K1))*   &
           DBLE(CF)**(1.0_dp/(1.0_dp+(LOG10(K1))**2))
   END FUNCTION FALL

   REAL(kind=dp) FUNCTION PAN ( A0,B0,C0,A1,B1,C1,CF)
      REAL A0,B0,C0,A1,B1,C1,CF
      REAL(kind=dp) K0, K1, KR, NC, F     
      K0 = DBLE(A0) * EXP(DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
      K1 = DBLE(A1) * EXP(DBLE(B1)/TEMP) * (TEMP/300.0_dp)**DBLE(C1)
      K0 = K0*CFACTOR
      KR = K0/K1
      NC = 0.75_dp-1.27_dp*(LOG10(DBLE(CF)))
      F  = 10.0_dp**(LOG10(DBLE(CF))/(1+(LOG10(KR)/NC)**2))
      PAN = K0*K1*F/(K0+K1)
   END FUNCTION PAN

   REAL(kind=dp) FUNCTION VILA ( A0,B0,C0,A1,B1,C1)
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) K0, K1     
      K0 = DBLE(A0) * EXP(DBLE(B0)/TEMP) * (TEMP/300.0_dp)**DBLE(C0)
      K1 = DBLE(A1) * EXP(DBLE(B1)/TEMP) * (TEMP/300.0_dp)**DBLE(C1)
      K0 = K0*CFACTOR
      VILA = K0*K1*0.35_dp/(K0+K1)
   END FUNCTION VILA

   REAL(kind=dp) FUNCTION NIT ( A0,B0,C0,n,X0,Y0)
      REAL A0,B0,C0,n,X0,Y0
      REAL(kind=dp) K0, K1, K2, K3, K4     
      K0 = 2.0E-22_dp * EXP(DBLE(n))
      K1 = 4.3E-1_dp*(TEMP/298.0_dp)**(-8)
      K0 = K0*CFACTOR
      K1 = K0/K1
      K2 = (K0/(1.0_dp+K1))*   &
           4.1E-1_dp**(1.0_dp/(1.0_dp+(LOG10(K1))**2))
      K3 = K2/(K2+DBLE(C0))
      K4 = DBLE(A0)*(DBLE(X0)-TEMP*DBLE(Y0))
      NIT = K4 * EXP(DBLE(B0)/TEMP) * K3
   END FUNCTION NIT

   REAL(kind=dp) FUNCTION ALK ( A0,B0,C0,n,X0,Y0)
      REAL A0,B0,C0,n,X0,Y0
      REAL(kind=dp) K0, K1, K2, K3, K4  
      K0 = 2.0E-22_dp * EXP(DBLE(n))
      K1 = 4.3E-1_dp*(TEMP/298.0_dp)**(-8)
      K0 = K0*CFACTOR
      K1 = K0/K1
      K2 = (K0/(1.0_dp+K1))*   &
           4.1E-1_dp**(1.0_dp/(1.0_dp+(LOG10(K1))**2))
      K3 = DBLE(C0)/(K2+DBLE(C0))
      K4 = DBLE(A0)*(DBLE(X0)-TEMP*DBLE(Y0))
      ALK = K4 * EXP(DBLE(B0)/TEMP) * K3
   END FUNCTION ALK

   REAL(kind=dp) FUNCTION EXTARR ( A0,B0,C0 )
      REAL A0,B0,C0      
      EXTARR =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (TEMP/300.0_dp)**DBLE(B0)
   END FUNCTION EXTARR

   REAL(kind=dp) FUNCTION IEMISS ( A0,B0,C0 )
      REAL A0,B0,C0 
      REAL(kind=dp) CT, CD, CL, BS, AT
      AT = DBLE(A0) * 8.314_dp * TEMP
      BS = DBLE(B0) * SUN
      CL = 2.88E-3_dp*BS/((1.0_dp+7.29E-6_dp*BS**2.0_dp)**0.5_dp)
      CT = EXP(9.5E+4_dp*(TEMP-DBLE(A0))/AT)
      CD = 1 + EXP(2.3E+5_dp*(TEMP-3.14E+2_dp)/AT)
      IEMISS = DBLE(C0) * CT * CL / CD
   END FUNCTION IEMISS

   REAL(kind=dp) FUNCTION MCMHV ( A0,B0,C0 )
      REAL A0,B0,C0
      IF (SUN.GT.0) THEN
         MCMHV = DBLE(A0) * SUN**DBLE(B0) * EXP(-DBLE(C0)/SUN)
      ELSE
         MCMHV = 0.0
      ENDIF
   END FUNCTION MCMHV

  !---------------------------------------------------------------------------

  ELEMENTAL REAL(kind=dp) FUNCTION k_3rd(temp,cair,k0_300K,n,kinf_300K,m,fc)

    INTRINSIC LOG10

    REAL(kind=dp), INTENT(IN) :: temp      ! temperature [K]
    REAL(kind=dp), INTENT(IN) :: cair      ! air concentration [molecules/cm3]
    REAL, INTENT(IN) :: k0_300K   ! low pressure limit at 300 K
    REAL, INTENT(IN) :: n         ! exponent for low pressure limit
    REAL, INTENT(IN) :: kinf_300K ! high pressure limit at 300 K
    REAL, INTENT(IN) :: m         ! exponent for high pressure limit
    REAL, INTENT(IN) :: fc        ! broadening factor (usually fc=0.6)
    REAL(kind=dp) :: zt_help, k0_T, kinf_T, k_ratio

    zt_help = 300._dp/temp
    k0_T    = k0_300K   * zt_help**(n) * cair ! k_0   at current T
    kinf_T  = kinf_300K * zt_help**(m)        ! k_inf at current T
    k_ratio = k0_T/kinf_T
    k_3rd   = k0_T/(1._dp+k_ratio)*fc**(1._dp/(1._dp+LOG10(k_ratio)**2))

  END FUNCTION k_3rd

  !---------------------------------------------------------------------------

  ELEMENTAL REAL(kind=dp) FUNCTION k_arr (k_298,tdep,temp)
    ! Arrhenius function

    REAL,     INTENT(IN) :: k_298 ! k at T = 298.15K
    REAL,     INTENT(IN) :: tdep  ! temperature dependence
    REAL(kind=dp), INTENT(IN) :: temp  ! temperature

    INTRINSIC EXP

    k_arr = k_298 * EXP(tdep*(1._dp/temp-3.3540E-3_dp)) ! 1/298.15=3.3540e-3

  END FUNCTION k_arr

  REAL(kind=dp) FUNCTION OH_O1D (J, H2O, TEMP, CFACTOR)
  REAL*8 J, H2O, TEMP, CFACTOR
  REAL*8 K1, K2, K3
  REAL*8 N2, O2

  N2=0.79d0
  O2=0.21d0

  K1=1.63d-10*exp(60d0/TEMP)
  K2=2.15d-11*exp(110d0/TEMP)
  K3=3.30d-11*exp(55d0/TEMP)

  OH_O1D=J*K1*H2O/(K1*H2O+K2*N2*CFACTOR+K3*O2*CFACTOR)
  END FUNCTION OH_O1D

  REAL(kind=dp) FUNCTION HO2_H2O(H2O,TEMP)
  REAL*8 TEMP, H2O
  HO2_H2O=1+1.4D-21*H2O*EXP(2200/TEMP)
  END FUNCTION HO2_H2O

  REAL(kind=dp) FUNCTION GCARR( A0,B0,C0 )
      REAL A0,B0,C0
      GCARR =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
  END FUNCTION GCARR

  REAL(kind=dp) FUNCTION GC_HO2NO3( A0,B0,C0,A1,B1,C1 )
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) :: R0,R1
      R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
      R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)

      GC_HO2NO3 = (R0+R1*CFACTOR)*(1.D0+1.4E-21_dp*CFACTOR*1.0E-2_dp &
                   * EXP(2200.E+0_dp/TEMP))
  END FUNCTION GC_HO2NO3

  REAL(kind=dp) FUNCTION GC_TBRANCH( A0,B0,C0,A1,B1,C1 )
! Temperature Dependent Branching Ratio                                                                                                     
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) :: R0,R1
      R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
      R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)

      GC_TBRANCH = R0/(1.D0+R1)

  END FUNCTION GC_TBRANCH

  REAL(kind=dp) FUNCTION GC_RO2HO2( A0,B0,C0,A1,B1,C1 )
! Carbon Dependence of RO2+HO2                                                                                                              
      REAL A0,B0,C0,A1,B1,C1
      REAL(kind=dp) :: R0,R1
      R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
      R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)

      GC_RO2HO2 = R0*(1E0-EXP(-0.245E0*R1))

  END FUNCTION GC_RO2HO2

  REAL(kind=dp) FUNCTION GC_DMSOH( A0,B0,C0,A1,B1,C1 )
! DMS+OH+O2                                                                                                                                 
    REAL A0,B0,C0,A1,B1,C1
    REAL(kind=dp) :: R0,R1
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)
!    GC_DMSOH = R0/(1e0_dp+R1*0.2095e0_dp)                                                                                                  
    GC_DMSOH = (R0*CFACTOR*0.2095e0_dp)/(1e0_dp+R1*0.2095e0_dp)

  END FUNCTION GC_DMSOH

  REAL(kind=dp) FUNCTION GC_GLYXNO3( A0,B0,C0 )
! ---  K = K1*([O2]+3.5D18)/(2*[O2]+3.5D18) --- HO2+2*CO branch of GLYX+OH/NO3                                                              
    REAL A0,B0,C0
    REAL(kind=dp) R0
    REAL(kind=dp) O2

    O2 = CFACTOR*0.2095e0_dp
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    GC_GLYXNO3 = R0*(O2+3.5E+18_dp)/(2.E+0_dp*O2+3.5E+18_dp)

  END FUNCTION GC_GLYXNO3

  REAL(kind=dp) FUNCTION GC_OHHNO3( A0,B0,C0,A1,B1,C1,A2,B2,C2 )
! ---  OH + HNO3:   K = K0 + K3[M] / (1 + K3[M]/K2)  ------                                                                                 
    REAL A0,B0,C0,A1,B1,C1,A2,B2,C2
    REAL(kind=dp) R0,R1,R2
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)
    R2 =  CFACTOR*(DBLE(A2) * EXP(DBLE(C2)/TEMP) * (300._dp/TEMP)**DBLE(B2))
    GC_OHHNO3 = R0 + R2/(1.E0_dp + R2/R1)

  END FUNCTION GC_OHHNO3

  REAL(kind=dp) FUNCTION GC_GLYCOHA( A0,B0,C0 )
!                                                                                                                                           
    REAL A0,B0,C0,R0,GLYC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    GLYC_FRAC=1e+0_dp-11.0729e+0_dp*EXP(-(1._dp/73._dp)*TEMP)
    IF (GLYC_FRAC<0e+0_dp) GLYC_FRAC=0e+0_dp
    GC_GLYCOHA=R0*GLYC_FRAC

  END FUNCTION GC_GLYCOHA

  REAL(kind=dp) FUNCTION GC_GLYCOHB( A0,B0,C0 )
!                                                                                                                                           
    REAL A0,B0,C0
    REAL(kind=dp) :: R0,GLYC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    GLYC_FRAC=1e+0_dp-11.0729e+0_dp*EXP(-(1._dp/73._dp)*TEMP)
    IF (GLYC_FRAC<0e+0_dp) GLYC_FRAC=0e+0_dp
    GC_GLYCOHB=R0*(1e0_dp-GLYC_FRAC)

  END FUNCTION GC_GLYCOHB

  REAL(kind=dp) FUNCTION GC_HACOHA( A0,B0,C0 )
!                                                                                                                                           
    REAL A0,B0,C0
    REAL(kind=dp) :: R0,HAC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    HAC_FRAC=1e+0_dp-23.7e+0_dp*EXP(-(1._dp/60._dp)*TEMP)
    IF (HAC_FRAC<0e+0_dp) HAC_FRAC=0e+0_dp
    GC_HACOHA=R0*HAC_FRAC

  END FUNCTION GC_HACOHA

  REAL(kind=dp) FUNCTION GC_HACOHB( A0,B0,C0 )
!                                                                                                                                           
    REAL A0,B0,C0
    REAL(kind=dp) :: R0,HAC_FRAC
    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    HAC_FRAC=1e+0_dp-23.7e+0_dp*EXP(-(1._dp/60._dp)*TEMP)
    IF (HAC_FRAC<0e+0_dp) HAC_FRAC=0e+0_dp
    GC_HACOHB=R0*(1.E0_dp-HAC_FRAC)

  END FUNCTION GC_HACOHB

  REAL(kind=dp) FUNCTION GC_OHCO( A0,B0,C0 )

    REAL A0,B0,C0,R0
    REAL KLO1,KLO2,KHI1,KHI2,XYRAT1,XYRAT2,BLOG1,BLOG2,FEXP1,FEXP2
    REAL KCO1,KCO2,KCO

    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R0 = R0 * (1.E+0_dp + 0.6e+0_dp*9.871E7_dp)

    ! new OH+CO rate from JPL2006.                                                                                                          
    KLO1=5.9E-33_dp*(300._dp/TEMP)**(1.E+0_dp)
    KHI1=1.1E-12_dp*(300._dp/TEMP)**(-1.3E0_dp)
    XYRAT1=KLO1*CFACTOR/KHI1
    BLOG1=LOG10(XYRAT1)
    FEXP1=1.E+0_dp/(1.E+0_dp+BLOG1*BLOG1)
    KCO1=KLO1*CFACTOR*0.6**FEXP1/(1.e+0_dp+XYRAT1)
    KLO2=1.5E-13_dp*(300._dp/TEMP)**(0.E+0_dp)
    KHI2=2.1e+09_dp *(300._dp/TEMP)**(-6.1E+0_dp)
    XYRAT2=KLO2*CFACTOR/KHI2
    BLOG2=LOG10(XYRAT2)
    FEXP2=1.E+0_dp/(1.E+0_dp+BLOG2*BLOG2)
    KCO2=KLO2*0.6**FEXP2/(1.e+0_dp+XYRAT2)
    KCO=KCO1+KCO2
    GC_OHCO=KCO

  END FUNCTION GC_OHCO

  REAL(kind=dp) FUNCTION GC_RO2NO( B,A0,B0,C0,A1,B1,C1 )
    ! ---  K = K1*(1-FYRNO3(K2,M,T))  ---  abstraction branch of RO2+NO                                                                     
    CHARACTER(*) B !Branch Toggle                                                                                                           
    REAL A0,B0,C0,A1,B1,C1
    REAL(kind=dp) :: R0,R1
    REAL(kind=dp) :: YYYN, XXYN,  AAA,  RARB, ZZYN
    REAL(kind=dp) :: XF, ALPHA, Y300, BETA, XMINF, XM0
    REAL(kind=dp) :: FYRNO3
    Y300 = 0.826
    ALPHA= 1.94e-22
    BETA = 0.97
    XM0  = 0.
    XMINF= 8.1
    XF   = 0.411

    R0 =  DBLE(A0) * EXP(DBLE(C0)/TEMP) * (300._dp/TEMP)**DBLE(B0)
    R1 =  DBLE(A1) * EXP(DBLE(C1)/TEMP) * (300._dp/TEMP)**DBLE(B1)

    ! Initialize static variables                                                                                                           

    XXYN   = ALPHA*EXP(BETA*R1)*CFACTOR*((300./TEMP)**XM0)
    YYYN   = Y300*((300./TEMP)**XMINF)
    AAA    = LOG10(XXYN/YYYN)
    ZZYN   = 1./(1.+ AAA*AAA )
    RARB   = (XXYN/(1.+ (XXYN/YYYN)))*(XF**ZZYN)
    FYRNO3 = RARB/(1. + RARB)
    IF (trim(B) .eq. 'A') THEN
       GC_RO2NO     = R0 * FYRNO3
    ELSEIF (trim(B) .eq. 'B') THEN
       GC_RO2NO     = R0 * (1.E+0_dp - FYRNO3)
    ENDIF

  END FUNCTION GC_RO2NO

  REAL(kind=dp) FUNCTION GCJPL3( k0_300, n, ki_300, m)
!  Functions given in JPL Booklet                                                                                                           
      REAL k0_300, n, ki_300,m
      REAL k0, ki

      k0=k0_300*((TEMP/300.d0)**(-n))
      ki=ki_300*((TEMP/300.d0)**(-m))

!      GCJPL3=(k0*CFACTOR)/(1+k0*CFACTOR/ki)*0.6** &                                                                                          
!       ((1+((LOG10(k0*CFACTOR/ki))**2d0)**-1.0d0))                                                                                          
      GCJPL3=(k0/(1.d0+k0/(ki/CFACTOR)))*0.6** &
        ((1+((log10(k0/(ki/CFACTOR)))**2d0)**1.0e0))
      GCJPL3=GCJPL3*CFACTOR
  END FUNCTION GCJPL3

  REAL(kind=dp) FUNCTION GCJPLEQ( A0,B0,C0,A1,B1,C1,A2,B2,C2,FV,FCT1,FCT2 )
 ! Function calculates the rate constant of the forward reaction                                                                            
 ! calculates the equilibrium constant                                                                                                      
 ! Find the backwards reaction by K=kforward/kbackwards                                                                                     
       REAL A0,B0,C0,A1,B1,C1
       REAL(kind=dp) :: R0,R1
       REAL, OPTIONAL :: A2,B2,C2,FV,FCT1,FCT2 !If a P-dependent rxn                                                                        

       ! Calculate Backwards reaction                                                                                                       
       R0 = GCARR( A0,B0,C0 )

       ! Calculate forwards reaction                                                                                                        
       IF (present(A2)) THEN ! P-dependent                                                                                                  
          IF (present(B2) .and. present(C2) .and. present(FV) &
                          .and. present(FCT1) .and. present(FCT2)) THEN
             R1 = GCJPLPR( A1,B1,C1,A2,B2,C2,FV,FCT1,FCT2)
          ELSE ! Missing params!                                                                                                            
             write(*,'(a)') 'GCJPLEQ: Missing parameters for P-dependent reaction.'
             write(*,'(a)') 'GCJPLEQ: Returning zero'
             GCJPLEQ = 0.E0
             RETURN
          ENDIF
       ELSE
          R1 = gcarr( A1,B1,C1 ) !Std. Arrhenius eqn.                                                                                       
       ENDIF

       GCJPLEQ=R1/R0
  END FUNCTION GCJPLEQ

  REAL(kind=dp) FUNCTION GCJPLPR(A0,B0,C0,A1,B1,C1,FV,FCT1,FCT2)
! * PRESSURE-DEPENDENT EFFECTS                                                                                                              
! * ADD THE THIRD BODY EFFECT FOR PRESSURE DEPENDENCE OF RATE                                                                               
! * COEFFICIENTS.                                                                                                                           
! A0 B0, & C0 are the Arrhenius parameters for the lower-limit                                                                              
! rate. A1, B1 & C1 are the upper-limit parameters.                                                                                         
! FV is the falloff curve paramter, (SEE ATKINSON ET. AL (1992)                                                                             
! J. PHYS. CHEM. REF. DATA 21, P. 1145). USUALLY = 0.6                                                                                      
!                                                                                                                                           
       REAL A0,B0,C0,A1,B1,C1,FV,FCT1,FCT2
       REAL FCT,XYRAT,BLOG,RLOW,RHIGH,FEXP

       RLOW  = GCARR( A0,B0,C0 )*CFACTOR
       RHIGH = GCARR( A1,B1,C1 )

       IF     (FCT2.NE.0.) THEN
             FCT            = EXP(-TEMP / FCT1) + EXP(-FCT2 / TEMP)
             XYRAT          = RLOW/RHIGH
             BLOG           = LOG10(XYRAT)
             FEXP           = 1.e+0_dp / (1.e+0_dp + BLOG * BLOG)
             GCJPLPR        = RLOW*FCT**FEXP/(1e+0_dp+XYRAT)
       ELSEIF (FCT1.NE.0.) THEN
             FCT            = EXP(-TEMP / FCT1)
             XYRAT          = RLOW/RHIGH
             BLOG           = LOG10(XYRAT)
             FEXP           = 1.e+0_dp / (1.e+0_dp + BLOG * BLOG)
             GCJPLPR        = RLOW*FCT**FEXP/(1e+0_dp+XYRAT)
       ELSE
             XYRAT          = RLOW/RHIGH
             BLOG           = LOG10(XYRAT)
             FEXP           = 1.e+0_dp / (1.e+0_dp + BLOG * BLOG)
             GCJPLPR        = RLOW*FV**FEXP/(1e+0_dp+XYRAT)
       ENDIF


  END FUNCTION GCJPLPR

  REAL(kind=dp) FUNCTION GCIUPAC3(ko_300,n,ki_300,m,Fc)
! Function calcualtes the rate constant of 3 body reaction using IUPAC                                                                      
! methology                                                                                                                                 
  REAL ko_300,n,ki_300,m,Fc
  REAL ko, ki, F, NN

  ko=ko_300*((TEMP/300.e0)**n)*CFACTOR
  ki=ki_300*((TEMP/300.e0)**m)

  NN=0.75-1.27*LOG10(Fc)
  F=10.0**(LOG10(Fc)/(1.0e0+(LOG10(ko/ki)/NN)**2.0))


  GCIUPAC3=ko/(1+ko/ki)*F
  END FUNCTION GCIUPAC3


!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!  End of User-defined Rate Law functions
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

! End Rate Law Functions from KPP_HOME/util/UserRateLaws


! Begin INLINED Rate Law Functions


! End INLINED Rate Law Functions

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! Update_SUN - update SUN light using TIME
!   Arguments :
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

  SUBROUTINE Update_SUN()
      !USE aromatics_kpp_Parameters
      !USE aromatics_kpp_Global

    IMPLICIT NONE

    REAL(kind=dp) :: SunRise, SunSet
    REAL(kind=dp) :: Thour, Tlocal, Ttmp 
    ! PI - Value of pi
    REAL(kind=dp), PARAMETER :: PI = 3.14159265358979d0
    
    SunRise = 4.5_dp 
    SunSet  = 19.5_dp 
    Thour = TIME/3600.0_dp 
    Tlocal = Thour - (INT(Thour)/24)*24

    IF ((Tlocal>=SunRise).AND.(Tlocal<=SunSet)) THEN
       Ttmp = (2.0*Tlocal-SunRise-SunSet)/(SunSet-SunRise)
       IF (Ttmp.GT.0) THEN
          Ttmp =  Ttmp*Ttmp
       ELSE
          Ttmp = -Ttmp*Ttmp
       END IF
       SUN = ( 1.0_dp + COS(PI*Ttmp) )/3.0_dp 
    ELSE
       SUN = 0.0_dp 
    END IF

 END SUBROUTINE Update_SUN

! End of Update_SUN function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! Update_RCONST - function to update rate constants
!   Arguments :
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

SUBROUTINE Update_RCONST ( )




! Begin INLINED RCONST

 

  RO2 = & 
  C(ind_OXYLO2) + C(ind_OXYBIPERO2) + C(ind_MXYLO2) + C(ind_MXYBIPERO2) + & 
   C(ind_PXYLO2) + C(ind_PXYBIPERO2) + C(ind_BZBIPERO2) + &  
   C(ind_C6H5CH2O2) + C(ind_TLBIPERO2) + C(ind_OXYMUCCO3) + &  
   C(ind_OXYMUCO2) + C(ind_MC6CO2O2) + C(ind_NOXYOLO2) + C(ind_OXYOLO2) + &  
   C(ind_MXYMUCCO3) + C(ind_MXYMUCO2) + C(ind_C726CO5O2) + C(ind_CH3CO3) &  
   + C(ind_MXYOLO2) + C(ind_NMXYOLO2) + C(ind_PXYMUCCO3) + &  
   C(ind_PXYMUCO2) + C(ind_C6M5CO2O2) + C(ind_NPXYOLO2) + C(ind_PXYOLO2) &  
   + C(ind_BZEMUCCO3) + C(ind_BZEMUCO2) + C(ind_C5DIALO2) + &  
   C(ind_NPHENO2) + C(ind_PHENO2) + C(ind_CRESO2) + C(ind_NCRESO2) + &  
   C(ind_TLEMUCCO3) + C(ind_TLEMUCO2) + C(ind_C615CO2O2) + C(ind_OXYLCO3) &  
   + C(ind_OXYL1O2) + C(ind_CH3O2) + C(ind_EPXM2DLCO3) + C(ind_C4MCO2O2) &  
   + C(ind_DM123O2) + C(ind_MXYLCO3) + C(ind_MXYL1O2) + C(ind_C3MCODBCO3) &  
   + C(ind_C5CO14O2) + C(ind_EPXMDLCO3) + C(ind_C3MDIALO2) + &  
   C(ind_MXY1O2) + C(ind_PXYLCO3) + C(ind_PXYL1O2) + C(ind_PXY1O2) + &  
   C(ind_MALDIALCO3) + C(ind_EPXDLCO3) + C(ind_C3DIALO2) + &  
   C(ind_MALDIALO2) + C(ind_C6H5O2) + C(ind_C6H5CO3) + C(ind_NBZFUO2) + &  
   C(ind_BZFUO2) + C(ind_BIACETO2) + C(ind_C5DICARBO2) + &  
   C(ind_MC4CODBCO3) + C(ind_MC5CO2OHO2) + C(ind_NOXYFUO2) + &  
   C(ind_C6OTKETO2) + C(ind_OXYFUO2) + C(ind_NTLFUO2) + C(ind_TLFUO2) + &  
   C(ind_OXCATEC1O2) + C(ind_C4MCODBCO3) + C(ind_C5MCO2OHO2) + &  
   C(ind_NMXYFUO2) + C(ind_C23O3MO2) + C(ind_MXYFUO2) + C(ind_MC3CODBCO3) &  
   + C(ind_C4M2ALOHO2) + C(ind_NPXYFUO2) + C(ind_MCOCOMOXO2) + &  
   C(ind_PXYFUO2) + C(ind_HCOCH2O2) + C(ind_MXCATEC1O2) + C(ind_DMKOHO2) &  
   + C(ind_C4CO2O2) + C(ind_PXCATEC1O2) + C(ind_HCOCOHCO3) + &  
   C(ind_CATEC1O2) + C(ind_MCATEC1O2) + C(ind_NTMB1FUO2) + &  
   C(ind_TMB1FUO2) + C(ind_HCOCO3) + C(ind_MECOACETO2) + C(ind_NOXYOL1O2) &  
   + C(ind_NOXYQO2) + C(ind_OXYQO2) + C(ind_OXNNCATCO2) + &  
   C(ind_OXNCATECO2) + C(ind_C534O2) + C(ind_CHOMOHCO3) + C(ind_CO2H3CO3) &  
   + C(ind_NMXYOL1O2) + C(ind_NMXYQO2) + C(ind_MXYQO2) + &  
   C(ind_MXNNCATCO2) + C(ind_MXNCATECO2) + C(ind_NPXYOL1O2) + &  
   C(ind_NPXYQO2) + C(ind_PXYQO2) + C(ind_PXNNCATCO2) + C(ind_PXNCATECO2) &  
   + C(ind_NPHEN1O2) + C(ind_NNCATECO2) + C(ind_NCATECO2) + C(ind_NBZQO2) &  
   + C(ind_PBZQO2) + C(ind_NPTLQO2) + C(ind_PTLQO2) + C(ind_NCRES1O2) + &  
   C(ind_MNNCATECO2) + C(ind_MNCATECO2) + C(ind_HOCH2CO3) + &  
   C(ind_CH3COCH2O2) + C(ind_MEKBO2) + C(ind_C6O4KETO2) + &  
   C(ind_NDNOXYOLO2) + C(ind_DNOXYOLO2) + C(ind_ACCOCOMEO2) + &  
   C(ind_C535O2) + C(ind_MMALANHYO2) + C(ind_CH3COCO3) + C(ind_MALANHYO2) &  
   + C(ind_NDNMXYOLO2) + C(ind_DNMXYOLO2) + C(ind_TL4OHNO2O2) + &  
   C(ind_NDNPXYOLO2) + C(ind_DNPXYOLO2) + C(ind_NDNPHENO2) + &  
   C(ind_DNPHENO2) + C(ind_NDNCRESO2) + C(ind_DNCRESO2) + &  
   C(ind_ACCOMECO3) + C(ind_C5M2OHOCO3) + C(ind_C4COMOHCO3) + &  
   C(ind_C23O3MCO3) + C(ind_C23O3CCO3) + C(ind_C7CO2OHCO3) + &  
   C(ind_C6MOHCOCO3) + C(ind_C5CO2OHCO3) + C(ind_C6CO2OHCO3) + &  
   C(ind_C4DBM2CO3) + C(ind_C5CO234O2) + C(ind_C5CO2DBCO3) + &  
   C(ind_C5DBCO2CO3) + C(ind_C4CO2DBCO3 ) + C(ind_MO2m) 


! End INLINED RCONST

  RCONST(1) = (7.45E4*(TEMP/300.)**-2.6)
  RCONST(2) = (8.0E-12*EXP(-2060./TEMP))
  RCONST(3) = (PAN(1.0E-31,0.0,-1.6,5.0E-11,0.0,-0.3,0.85))
  RCONST(4) = (5.5E-12*EXP(188./TEMP))
  RCONST(5) = (PAN(1.3E-31,0.0,-1.5,2.3E-11,0.0,0.24,0.6))
  RCONST(6) = (1.68E8*EXP(67./TEMP)+3.95E8*EXP(130./TEMP))
  RCONST(7) = (1.4E-12*EXP(-1310./TEMP))
  RCONST(8) = (1.4E-13*EXP(-2470./TEMP))
  RCONST(9) = (1.7E-20*EXP(530./TEMP))
  RCONST(10) = (1.8E-11*EXP(110./TEMP))
  RCONST(11) = (4.50E-14*EXP(-1260./TEMP))
  RCONST(12) = (PAN(3.6E-30,0.0,-4.1,1.9E-12,0.0,0.2,0.35))
! RCONST(13) = constant rate coefficient
  RCONST(14) = (1.70E-12*EXP(-940./TEMP))
  RCONST(15) = (7.7E-12*EXP(-2100./TEMP))
  RCONST(16) = (2.9E-12*EXP(-160./TEMP))
  RCONST(17) = (2.03E-16*EXP(693./TEMP)*(TEMP/300.)**4.57)
  RCONST(18) = (4.8E-11*EXP(250./TEMP))
  RCONST(19) = ((1+(3.5E-4*EXP(2200./TEMP)))*(2.20E-13*EXP(600./TEMP)+4.75E-14*EXP(980./TEMP)))
  RCONST(20) = (PAN(7.4E-31,0.0,-2.4,3.3E-11,0.0,-0.3,0.81))
  RCONST(21) = (PAN(3.2E-30,0.0,-4.5,3.0E-11,0.0,0.0,0.41))
! RCONST(22) = constant rate coefficient
  RCONST(23) = (3.45E-12*EXP(270./TEMP))
  RCONST(24) = (PAN(1.4E-31,0.0,-3.1,4.0E-12,0.0,0.0,0.4))
  RCONST(25) = (3.2E-13*EXP(690./TEMP))
! RCONST(26) = constant rate coefficient
  RCONST(27) = (2.5E-12*EXP(260./TEMP))
  RCONST(28) = (2.4E-14*EXP(460./TEMP)+(1.625E-14*EXP(1335./TEMP))/(1+(602.*EXP(-864./TEMP))))
  RCONST(29) = (SUN*3.78E-5)
  RCONST(30) = (SUN*4.41E-4)
  RCONST(31) = (SUN*7.88E-6)
  RCONST(32) = (SUN*8.92E-3)
  RCONST(33) = (SUN*2.23E-2)
  RCONST(34) = (SUN*1.54E-1)
  RCONST(35) = (SUN*1.98E-3)
  RCONST(36) = (SUN*6.85E-7)
  RCONST(37) = (GCARR(8.50E-13,0.0E+00,-2450.0))
  RCONST(38) = (SUN*4.5E-5)
  RCONST(39) = (SUN*5.21E-6)
  RCONST(40) = (PAN(1.3E-3,-1.1E4,-3.5,9.7E14,-1.108E4,0.1,0.35))
  RCONST(41) = (PAN(4.1E-5,-10650.0,0.0,6.0E15,-11170.0,0.0,0.4))
  RCONST(42) = (GCARR(2.80E-12,0.0E+00,300.0))
  RCONST(43) = (GCARR(4.10E-13,0.0E+00,750.0))
  RCONST(44) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,2.62E1,0.0,-1130.0))
  RCONST(45) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,4.E-2,0.0,1130.0))
  RCONST(46) = (GCJPLPR(1.00E-30,4.8E+00,0.0,7.2E-12,2.1E0,0.0,0.6,0.0,0.0))
  RCONST(47) = (GCJPLPR(1.05E-02,4.8E+00,-11234.0,7.58E16,2.1E0,-11234.0,0.6,0.0,0.0))
  RCONST(48) = (GCARR(2.66E-12,0.0E+00,200.0))
  RCONST(49) = (GCARR(1.14E-12,0.0E+00,200.0))
  RCONST(50) = (GCARR(5.50E-12,0.0E+00,125.0))
  RCONST(51) = (GCARR(5.80E-16,0.0E+00,0.0))
  RCONST(52) = (GCARR(2.90E-12,0.0E+00,-345.0))
  RCONST(53) = (KCO(1.44E-13,4.2E19))
  RCONST(54) = (SUN*6.49E-6)
  RCONST(55) = (SUN*4.65E-5)
  RCONST(56) = (SUN*6.01E-5)
  RCONST(57) = (SUN*0.05*1.95E-4)
  RCONST(58) = (SUN*0.95*1.95E-4)
  RCONST(59) = (GCARR(2.45E-12,0.0E+00,-1775.0))
  RCONST(60) = (GCARR(2.80E-12,0.0E+00,300.0))
  RCONST(61) = (GCARR(4.10E-13,0.0E+00,750.0))
  RCONST(62) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,2.62E1,0.0,-1130.0))
  RCONST(63) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,4.E-2,0.0,1130.0))
  RCONST(64) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,2.62E1,0.0,-1130.0))
  RCONST(65) = (GC_TBRANCH(9.50E-14,0.0E+00,390.0,4.E-2,0.0,1130.0))
  RCONST(66) = (GCJPLPR(1.00E-30,4.8E+00,0.0,7.2E-12,2.1E0,0.0,0.6,0.0,0.0))
  RCONST(67) = (GCJPLPR(1.05E-02,4.8E+00,-11234.0,7.58E16,2.1E0,-11234.0,0.6,0.0,0.0))
  RCONST(68) = (GCARR(2.66E-12,0.0E+00,200.0))
  RCONST(69) = (GCARR(1.14E-12,0.0E+00,200.0))
  RCONST(70) = (GCARR(5.50E-12,0.0E+00,125.0))
  RCONST(71) = (GCARR(5.80E-16,0.0E+00,0.0))
  RCONST(72) = (GCARR(2.90E-12,0.0E+00,-345.0))
  RCONST(73) = (KCO(1.44E-13,4.2E19))
  RCONST(74) = (SUN*6.49E-6)
  RCONST(75) = (SUN*4.65E-5)
  RCONST(76) = (SUN*6.01E-5)
  RCONST(77) = (SUN*0.05*1.95E-4)
  RCONST(78) = (SUN*0.95*1.95E-4)
  RCONST(79) = (SUN*3.1E-3)
  RCONST(80) = (GCARR(3.36E-12,0.0E+00,-1860.0))
  RCONST(81) = (GCARR(1.50E-11,0.0E+00,0.0))
  RCONST(82) = (SUN*8.46E-5)
  RCONST(83) = (SUN*8.61E-6)
  RCONST(84) = (SUN*3.71E-5)
  RCONST(85) = (GC_GLYXNO3(1.40E-12,0.0E+00,-1860.0))
  RCONST(86) = (GCARR(3.10E-12,0.0E+00,340.0))
  RCONST(87) = (SUN*6.4E-6)
  RCONST(88) = (GCARR(4.63E-12,0.0E+00,350.0))
  RCONST(89) = (GCARR(1.40E-12,0.0E+00,-1900.0))
  RCONST(90) = (GCARR(6.13E-13,0.0E+00,200.0))
  RCONST(91) = (SUN*3.0E-5)
  RCONST(92) = (GCJPLEQ(9.30E-29,0.0E+00,14000.0,9.7E-29,5.6E0,0.0,9.3E-12,1.5E0,0.,0.6,0.,0.))
  RCONST(93) = (SUN*1.0E-6)
  RCONST(94) = (GCARR(1.80E-12,0.0E+00,500.0)*0.1*RO2)
  RCONST(95) = (GCARR(1.80E-12,0.0E+00,500.0)*0.9*RO2)
  RCONST(96) = (GCJPLPR(9.70E-29,5.6E+00,0.0,9.3E-12,1.5E0,0.0,0.6,0.0,0.0))
  RCONST(97) = (GCARR(8.10E-12,0.0E+00,270.0))
  RCONST(98) = (GCARR(5.20E-13,0.0E+00,980.0))
  RCONST(99) = (GC_GLYCOHA(8.00E-12,0.0E+00,0.0))
  RCONST(100) = (GC_GLYCOHB(8.00E-12,0.0E+00,0.0))
  RCONST(101) = (SUN*1.67E-5)
  RCONST(102) = (SUN*1.2e-5)
  RCONST(103) = (SUN*1.2e-5)
  RCONST(104) = (SUN*1.2e-5)
  RCONST(105) = (SUN*1.2e-5)
  RCONST(106) = (SUN*1.2e-5)
  RCONST(107) = (SUN*1.2e-5)
  RCONST(108) = (SUN*1.2e-5)
  RCONST(109) = (SUN*1.2e-5)
  RCONST(110) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(111) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(112) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(113) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(114) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(115) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(116) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(117) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(118) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(119) = (k_3rd(TEMP,2.5E+19_dp,9.7E-29,5.6,9.3E-12,1.5,0.6))
  RCONST(120) = (1.1D-11*EXP(364/TEMP)*0.65)
  RCONST(121) = (1.1D-11*EXP(364/TEMP)*0.20)
  RCONST(122) = (1.1D-11*EXP(364/TEMP)*0.15)
  RCONST(123) = (2.91D-13*EXP(1300/TEMP)*0.77*0.6)
  RCONST(124) = (2.91D-13*EXP(1300/TEMP)*0.77*0.4)
  RCONST(125) = (2.91D-13*EXP(1300/TEMP)*0.82*0.6)
  RCONST(126) = (2.91D-13*EXP(1300/TEMP)*0.82*0.4)
  RCONST(127) = (2.91D-13*EXP(1300/TEMP)*0.859*0.6)
  RCONST(128) = (2.91D-13*EXP(1300/TEMP)*0.859*0.4)
  RCONST(129) = (2.91D-13*EXP(1300/TEMP)*0.859*0.6)
  RCONST(130) = (2.91D-13*EXP(1300/TEMP)*0.859*0.4)
  RCONST(131) = (2.91D-13*EXP(1300/TEMP)*0.859*0.6)
  RCONST(132) = (2.91D-13*EXP(1300/TEMP)*0.859*0.4)
! RCONST(133) = constant rate coefficient
! RCONST(134) = constant rate coefficient
! RCONST(135) = constant rate coefficient
! RCONST(136) = constant rate coefficient
! RCONST(137) = constant rate coefficient
  RCONST(138) = (2.3D-12*EXP(-190/TEMP)*0.352)
  RCONST(139) = (2.3D-12*EXP(-190/TEMP)*0.118)
  RCONST(140) = (2.3D-12*EXP(-190/TEMP)*0.53)
  RCONST(141) = (1.8D-12*EXP(340/TEMP)*0.07)
  RCONST(142) = (1.8D-12*EXP(340/TEMP)*0.18)
  RCONST(143) = (1.8D-12*EXP(340/TEMP)*0.65)
  RCONST(144) = (1.8D-12*EXP(340/TEMP)*0.10)
  RCONST(145) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(146) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(147) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(148) = constant rate coefficient
  RCONST(149) = (1.30D-12*0.20*RO2)
  RCONST(150) = (1.30D-12*0.60*RO2)
  RCONST(151) = (1.30D-12*0.20*RO2)
  RCONST(152) = (2.7D-12*EXP(360/TEMP)*0.138)
  RCONST(153) = (2.7D-12*EXP(360/TEMP)*0.862)
! RCONST(154) = constant rate coefficient
  RCONST(155) = (9.20D-14*RO2*0.70)
  RCONST(156) = (9.20D-14*RO2*0.30)
  RCONST(157) = (2*1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(158) = (2.00D-18)
  RCONST(159) = (5.99D-11*0.632)
  RCONST(160) = (5.99D-11*0.368)
  RCONST(161) = (SUN*8.92E-3*0.05*0.5)
  RCONST(162) = (SUN*8.92E-3*0.05*0.5)
  RCONST(163) = (3.20D-11*0.51)
  RCONST(164) = (3.20D-11*0.39)
  RCONST(165) = (3.20D-11*0.10)
  RCONST(166) = (8.00D-11*0.07)
  RCONST(167) = (8.00D-11*0.73)
  RCONST(168) = (8.00D-11*0.20)
  RCONST(169) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(170) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(171) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(172) = constant rate coefficient
  RCONST(173) = (1.30D-12*0.20*RO2)
  RCONST(174) = (1.30D-12*0.60*RO2)
  RCONST(175) = (1.30D-12*0.20*RO2)
  RCONST(176) = (2.7D-12*EXP(360/TEMP)*0.138)
  RCONST(177) = (2.7D-12*EXP(360/TEMP)*0.862)
! RCONST(178) = constant rate coefficient
  RCONST(179) = (8.80D-13*RO2*0.6)
  RCONST(180) = (8.80D-13*RO2*0.2)
  RCONST(181) = (8.80D-13*RO2*0.2)
  RCONST(182) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(183) = (2.00D-18)
  RCONST(184) = (8.02D-11*0.26)
  RCONST(185) = (8.02D-11*0.74)
  RCONST(186) = (SUN*8.92E-3*0.028*0.5)
  RCONST(187) = (SUN*8.92E-3*0.028*0.5)
  RCONST(188) = (4.50D-11*0.39)
  RCONST(189) = (4.50D-11*0.10)
  RCONST(190) = (4.50D-11*0.51)
  RCONST(191) = (9.10D-11*0.07)
  RCONST(192) = (9.10D-11*0.42)
  RCONST(193) = (9.10D-11*0.51)
  RCONST(194) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(195) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(196) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(197) = constant rate coefficient
  RCONST(198) = (1.30D-12*0.20*RO2)
  RCONST(199) = (1.30D-12*0.60*RO2)
  RCONST(200) = (1.30D-12*0.20*RO2)
  RCONST(201) = (2.7D-12*EXP(360/TEMP)*0.138)
  RCONST(202) = (2.7D-12*EXP(360/TEMP)*0.862)
! RCONST(203) = constant rate coefficient
  RCONST(204) = (9.20D-14*RO2*0.70)
  RCONST(205) = (9.20D-14*RO2*0.30)
  RCONST(206) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(207) = (5.00D-18)
  RCONST(208) = (7.88D-11*0.26)
  RCONST(209) = (7.88D-11*0.74)
  RCONST(210) = (SUN*8.92E-3*0.1*0.5)
  RCONST(211) = (SUN*8.92E-3*0.1*0.5)
  RCONST(212) = (3.48D-11*0.51)
  RCONST(213) = (3.48D-11*0.39)
  RCONST(214) = (3.48D-11*0.10)
  RCONST(215) = (8.00D-11*0.07)
  RCONST(216) = (8.00D-11*0.73)
  RCONST(217) = (8.00D-11*0.2)
  RCONST(218) = (2.7D-12*EXP(360/TEMP)*0.082)
  RCONST(219) = (2.7D-12*EXP(360/TEMP)*0.918)
! RCONST(220) = constant rate coefficient
  RCONST(221) = (8.80D-13*RO2*0.2)
  RCONST(222) = (8.80D-13*RO2*0.6)
  RCONST(223) = (8.80D-13*RO2*0.2)
  RCONST(224) = (2*1.4D-12*EXP(-1860/TEMP)*2.75)
  RCONST(225) = (2.00D-18)
  RCONST(226) = (6.08D-11*0.31)
  RCONST(227) = (6.08D-11*0.69)
  RCONST(228) = (SUN*8.92E-3*0.1*0.5)
  RCONST(229) = (SUN*8.92E-3*0.1*0.5)
  RCONST(230) = (3.8D-12*0.742)
  RCONST(231) = (3.8D-12*0.258)
  RCONST(232) = (4.7D-13*EXP(1220/TEMP)*0.06)
  RCONST(233) = (4.7D-13*EXP(1220/TEMP)*0.8)
  RCONST(234) = (4.7D-13*EXP(1220/TEMP)*0.14)
  RCONST(235) = (1.5D-13*EXP(1310/TEMP))
  RCONST(236) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(237) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(238) = constant rate coefficient
  RCONST(239) = (2*(1.03D-13*EXP(365/TEMP)*2.4D-14*EXP(1620/TEMP))**0.5*RO2*0.2)
  RCONST(240) = (2*(1.03D-13*EXP(365/TEMP)*2.4D-14*EXP(1620/TEMP))**0.5*RO2*0.6)
  RCONST(241) = (2*(1.03D-13*EXP(365/TEMP)*2.4D-14*EXP(1620/TEMP))**0.5*RO2*0.2)
  RCONST(242) = (1.4D-11*0.103)
  RCONST(243) = (1.4D-11*0.506)
  RCONST(244) = (1.4D-11*0.391)
  RCONST(245) = (4.65D-11*0.200)
  RCONST(246) = (4.65D-11*0.727)
  RCONST(247) = (4.65D-11*0.073)
  RCONST(248) = (2.7D-12*EXP(360/TEMP)*0.889)
  RCONST(249) = (2.7D-12*EXP(360/TEMP)*0.111)
! RCONST(250) = constant rate coefficient
  RCONST(251) = (8.80D-13*RO2*0.20)
  RCONST(252) = (8.80D-13*RO2*0.60)
  RCONST(253) = (8.80D-13*RO2*0.20)
  RCONST(254) = (1.4D-12*EXP(-1860/TEMP)*2.75)
  RCONST(255) = (5.00D-18)
  RCONST(256) = (7.99D-11*0.31)
  RCONST(257) = (7.99D-11*0.69)
  RCONST(258) = (SUN*8.92E-3*0.1*0.5)
  RCONST(259) = (SUN*8.92E-3*0.1*0.5)
  RCONST(260) = (2.75D-11)
  RCONST(261) = (SUN*5.79E-6)
  RCONST(262) = (1.31D-11)
  RCONST(263) = (SUN*1.79E-6)
  RCONST(264) = (2.50D-14*EXP(-300/TEMP))
  RCONST(265) = (1.4D-12*EXP(-1860/TEMP))
  RCONST(266) = (1.29D-11)
  RCONST(267) = (SUN*1.10E-6)
  RCONST(268) = (SUN*1.10E-6)
  RCONST(269) = (1.60D-11)
  RCONST(270) = (7.62D-11)
  RCONST(271) = (SUN*5.79E-6)
  RCONST(272) = (7.28D-11)
  RCONST(273) = (SUN*8.33E-6)
  RCONST(274) = (1.00D+06*0.175)
  RCONST(275) = (1.00D+06*0.245)
  RCONST(276) = (1.00D+06*0.175)
  RCONST(277) = (1.00D+06*0.08)
  RCONST(278) = (1.00D+06*0.08)
  RCONST(279) = (1.00D+06*0.245)
  RCONST(280) = (1.11D-10)
  RCONST(281) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(282) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(283) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(284) = (7.5D-12*EXP(290/TEMP))
  RCONST(285) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(286) = (2.3E-12*1.74)
  RCONST(287) = (1.00D-11*RO2*0.70)
  RCONST(288) = (1.00D-11*RO2*0.30)
  RCONST(289) = (2*1.4D-12*EXP(-1860/TEMP)*8.5)
  RCONST(290) = (4.24D-11)
  RCONST(291) = (SUN*5.50E-5*2)
  RCONST(292) = (1.00D+06*0.18)
  RCONST(293) = (1.00D+06*0.125)
  RCONST(294) = (1.00D+06*0.125)
  RCONST(295) = (1.00D+06*0.57)
  RCONST(296) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(297) = (2.7D-12*EXP(360/TEMP)*0.173)
  RCONST(298) = (2.7D-12*EXP(360/TEMP)*0.827)
! RCONST(299) = constant rate coefficient
  RCONST(300) = (8.80D-12*RO2*0.20)
  RCONST(301) = (8.80D-12*RO2*0.60)
  RCONST(302) = (8.80D-12*RO2*0.20)
  RCONST(303) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(304) = (2.7D-12*EXP(360/TEMP))
! RCONST(305) = constant rate coefficient
  RCONST(306) = (9.20D-14*RO2*0.70)
  RCONST(307) = (9.20D-14*RO2*0.30)
  RCONST(308) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(309) = (2.7D-12*EXP(360/TEMP))
! RCONST(310) = constant rate coefficient
  RCONST(311) = (8.00D-13*RO2*0.70)
  RCONST(312) = (8.00D-13*RO2*0.30)
  RCONST(313) = (2.08D-12)
  RCONST(314) = (2.86D-13)
  RCONST(315) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(316) = (2.7D-12*EXP(360/TEMP))
! RCONST(317) = constant rate coefficient
  RCONST(318) = (8.00D-13*RO2*0.70)
  RCONST(319) = (8.00D-13*RO2*0.30)
  RCONST(320) = (2.01D-10)
  RCONST(321) = (2.81D-17)
  RCONST(322) = (2.05D-10)
  RCONST(323) = (2.05D-11)
  RCONST(324) = (SUN*5.79E-6)
  RCONST(325) = (6.03D-12)
  RCONST(326) = (SUN*1.79E-6)
  RCONST(327) = (2.50D-14*EXP(-300/TEMP))
  RCONST(328) = (1.4D-12*EXP(-1860/TEMP))
  RCONST(329) = (1.29D-11)
  RCONST(330) = (SUN*1.10E-6)
  RCONST(331) = (SUN*1.10E-6)
  RCONST(332) = (8.93D-12)
  RCONST(333) = (1.17D-10)
  RCONST(334) = (SUN*5.79E-6)
  RCONST(335) = (7.15D-11)
  RCONST(336) = (SUN*2.99E-6)
  RCONST(337) = (1.00D+06*0.11)
  RCONST(338) = (1.00D+06*0.11)
  RCONST(339) = (1.00D+06*0.26)
  RCONST(340) = (1.00D+06*0.26)
  RCONST(341) = (1.00D+06*0.26)
  RCONST(342) = (2.50D-10)
  RCONST(343) = (7.93D-11*0.5)
  RCONST(344) = (7.93D-11*0.5)
  RCONST(345) = (SUN*3.98E-6*0.5)
  RCONST(346) = (SUN*3.98E-6*0.5)
  RCONST(347) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(348) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(349) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(350) = (7.5D-12*EXP(290/TEMP))
  RCONST(351) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(352) = (2.3E-12*1.74)
  RCONST(353) = (1.00D-11*RO2*0.70)
  RCONST(354) = (1.00D-11*RO2*0.30)
  RCONST(355) = (2*1.4D-12*EXP(-1860/TEMP)*8.5)
  RCONST(356) = (4.28D-11)
  RCONST(357) = (SUN*5.50E-5*2)
  RCONST(358) = (1.00D+06*0.125)
  RCONST(359) = (1.00D+06*0.57)
  RCONST(360) = (1.00D+06*0.125)
  RCONST(361) = (1.00D+06*0.18)
  RCONST(362) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(363) = (2.7D-12*EXP(360/TEMP)*0.173)
  RCONST(364) = (2.7D-12*EXP(360/TEMP)*0.827)
! RCONST(365) = constant rate coefficient
  RCONST(366) = (8.80D-13*RO2*0.2)
  RCONST(367) = (8.80D-13*RO2*0.6)
  RCONST(368) = (8.80D-13*RO2*0.2)
  RCONST(369) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(370) = (2.7D-12*EXP(360/TEMP))
! RCONST(371) = constant rate coefficient
  RCONST(372) = (8.80D-13*RO2*0.2)
  RCONST(373) = (8.80D-13*RO2*0.2)
  RCONST(374) = (8.80D-13*RO2*0.6)
  RCONST(375) = (2.08D-12)
  RCONST(376) = (2.86D-13)
  RCONST(377) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(378) = (2.7D-12*EXP(360/TEMP))
! RCONST(379) = constant rate coefficient
  RCONST(380) = (8.00D-13*RO2*0.70)
  RCONST(381) = (8.00D-13*RO2*0.30)
  RCONST(382) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(383) = (2.7D-12*EXP(360/TEMP))
! RCONST(384) = constant rate coefficient
  RCONST(385) = (8.00D-13*RO2*0.70)
  RCONST(386) = (8.00D-13*RO2*0.30)
  RCONST(387) = (2.01D-10)
  RCONST(388) = (2.81D-17)
  RCONST(389) = (2.05D-10)
  RCONST(390) = (2.05D-11)
  RCONST(391) = (SUN*5.79E-6)
  RCONST(392) = (6.03D-12)
  RCONST(393) = (SUN*1.79E-6)
  RCONST(394) = (2.50D-14*EXP(-300/TEMP))
  RCONST(395) = (1.4D-12*EXP(-1860/TEMP))
  RCONST(396) = (1.29D-11)
  RCONST(397) = (SUN*1.10E-6)
  RCONST(398) = (SUN*1.10E-6)
  RCONST(399) = (8.93D-12)
  RCONST(400) = (7.65D-11)
  RCONST(401) = (SUN*5.79E-6)
  RCONST(402) = (7.16D-11)
  RCONST(403) = (SUN*8.33E-6)
  RCONST(404) = (1.00D+06*0.64)
  RCONST(405) = (1.00D+06*0.18)
  RCONST(406) = (1.00D+06*0.18)
  RCONST(407) = (1.12D-10)
  RCONST(408) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(409) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(410) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(411) = (7.5D-12*EXP(290/TEMP))
  RCONST(412) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(413) = (2.3E-12*1.74)
  RCONST(414) = (1.00D-11*RO2*0.70)
  RCONST(415) = (1.00D-11*RO2*0.30)
  RCONST(416) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(417) = (2.7D-12*EXP(360/TEMP)*0.173)
  RCONST(418) = (2.7D-12*EXP(360/TEMP)*0.827)
! RCONST(419) = constant rate coefficient
  RCONST(420) = (8.80D-12*RO2*0.2)
  RCONST(421) = (8.80D-12*RO2*0.6)
  RCONST(422) = (8.80D-12*RO2*0.2)
  RCONST(423) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(424) = (2.7D-12*EXP(360/TEMP))
! RCONST(425) = constant rate coefficient
  RCONST(426) = (9.20D-14*RO2*0.70)
  RCONST(427) = (9.20D-14*RO2*0.30)
  RCONST(428) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(429) = (2.7D-12*EXP(360/TEMP))
! RCONST(430) = constant rate coefficient
  RCONST(431) = (8.00D-13*RO2*0.70)
  RCONST(432) = (8.00D-13*RO2*0.30)
  RCONST(433) = (2.08D-12)
  RCONST(434) = (2.86D-13)
  RCONST(435) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(436) = (2.7D-12*EXP(360/TEMP))
! RCONST(437) = constant rate coefficient
  RCONST(438) = (8.00D-13*RO2*0.70)
  RCONST(439) = (8.00D-13*RO2*0.30)
  RCONST(440) = (2.01D-10)
  RCONST(441) = (2.81D-17)
  RCONST(442) = (2.05D-10)
  RCONST(443) = (9.77D-11)
  RCONST(444) = (SUN*5.79E-6)
  RCONST(445) = (7.30D-11)
  RCONST(446) = (SUN*2.99E-6)
  RCONST(447) = (1.00D+06*0.50)
  RCONST(448) = (1.00D+06*0.50)
  RCONST(449) = (1.21D-10)
  RCONST(450) = (8.16D-11)
  RCONST(451) = (SUN*3.98E-6)
  RCONST(452) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(453) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(454) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(455) = (7.5D-12*EXP(290/TEMP))
  RCONST(456) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(457) = (2.3E-12*1.74)
  RCONST(458) = (1.00D-11*RO2*0.3)
  RCONST(459) = (1.00D-11*RO2*0.7)
  RCONST(460) = (2*1.4D-12*EXP(-1860/TEMP)*4.0)
  RCONST(461) = (4.32D-11)
  RCONST(462) = (SUN*5.50E-5*2)
  RCONST(463) = (2.91D-13*EXP(1300/TEMP)*0.77)
  RCONST(464) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(465) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(466) = constant rate coefficient
  RCONST(467) = (8.80D-13*0.2*RO2)
  RCONST(468) = (8.80D-13*0.6*RO2)
  RCONST(469) = (8.80D-13*0.2*RO2)
  RCONST(470) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(471) = (2.7D-12*EXP(360/TEMP))
! RCONST(472) = constant rate coefficient
  RCONST(473) = (8.80D-13*RO2*0.2)
  RCONST(474) = (8.80D-13*RO2*0.6)
  RCONST(475) = (8.80D-13*RO2*0.2)
  RCONST(476) = (2*1.4D-12*EXP(-1860/TEMP)*2.0)
  RCONST(477) = (2.00D-18)
  RCONST(478) = (5.20D-11*0.83)
  RCONST(479) = (5.20D-11*0.17)
  RCONST(480) = (SUN*8.92E-3*0.14*0.4)
  RCONST(481) = (SUN*8.92E-3*0.14*0.6)
  RCONST(482) = (2.08D-12)
  RCONST(483) = (2.86D-13)
  RCONST(484) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(485) = (2.7D-12*EXP(360/TEMP))
! RCONST(486) = constant rate coefficient
  RCONST(487) = (8.00D-13*RO2*0.7)
  RCONST(488) = (8.00D-13*RO2*0.3)
  RCONST(489) = (9.9D-11)
  RCONST(490) = (9.2D-18)
  RCONST(491) = (1.0D-10)
  RCONST(492) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(493) = (2.7D-12*EXP(360/TEMP))
! RCONST(494) = constant rate coefficient
  RCONST(495) = (8.00D-13*RO2*0.7)
  RCONST(496) = (8.00D-13*RO2*0.3)
  RCONST(497) = (2.05D-11)
  RCONST(498) = (SUN*5.79E-6)
  RCONST(499) = (6.03D-12)
  RCONST(500) = (SUN*1.79E-6)
  RCONST(501) = (2.50D-14*EXP(-300/TEMP))
  RCONST(502) = (SUN*3.30E-6)
  RCONST(503) = (2.40D-15)
  RCONST(504) = (5.9D-12*EXP(225/TEMP))
  RCONST(505) = (2.7D-11)
  RCONST(506) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(507) = (2.7D-12*EXP(360/TEMP))
! RCONST(508) = constant rate coefficient
  RCONST(509) = (8.00D-13*RO2*0.70)
  RCONST(510) = (8.00D-13*RO2*0.30)
  RCONST(511) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(512) = (2.7D-12*EXP(360/TEMP))
! RCONST(513) = constant rate coefficient
  RCONST(514) = (8.00D-13*RO2*0.70)
  RCONST(515) = (8.00D-13*RO2*0.30)
  RCONST(516) = (2.08D-12)
  RCONST(517) = (2.86D-13)
  RCONST(518) = (1.7D-10*1.0)
  RCONST(519) = (2.8D-17)
  RCONST(520) = (2.0D-10*1.0)
  RCONST(521) = (9.64D-11)
  RCONST(522) = (SUN*5.79E-6)
  RCONST(523) = (1.00D+06*0.2)
  RCONST(524) = (1.00D+06*0.2)
  RCONST(525) = (1.00D+06*0.2)
  RCONST(526) = (1.00D+06*0.2)
  RCONST(527) = (1.00D+06*0.2)
  RCONST(528) = (7.16D-11)
  RCONST(529) = (SUN*2.99E-6)
  RCONST(530) = (1.20D-10)
  RCONST(531) = (7.99D-11)
  RCONST(532) = (SUN*3.98E-6)
  RCONST(533) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(534) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(535) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(536) = (7.5D-12*EXP(290/TEMP))
  RCONST(537) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(538) = (2.3E-12*1.74)
  RCONST(539) = (1.00D-11*RO2*0.7)
  RCONST(540) = (1.00D-11*RO2*0.3)
  RCONST(541) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(542) = (2.7D-12*EXP(360/TEMP)*0.105)
  RCONST(543) = (2.7D-12*EXP(360/TEMP)*0.895)
! RCONST(544) = constant rate coefficient
  RCONST(545) = (8.80D-13*RO2*0.20)
  RCONST(546) = (8.80D-13*RO2*0.60)
  RCONST(547) = (8.80D-13*RO2*0.20)
  RCONST(548) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(549) = (2.7D-12*EXP(360/TEMP))
! RCONST(550) = constant rate coefficient
  RCONST(551) = (8.80D-13*RO2*0.2)
  RCONST(552) = (8.80D-13*RO2*0.6)
  RCONST(553) = (8.80D-13*RO2*0.2)
  RCONST(554) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(555) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(556) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(557) = (7.5D-12*EXP(290/TEMP))
  RCONST(558) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(559) = (2.3E-12*1.74)
  RCONST(560) = (1.00D-11*0.70*RO2)
  RCONST(561) = (1.00D-11*0.30*RO2)
  RCONST(562) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(563) = (2.7D-12*EXP(360/TEMP))
! RCONST(564) = constant rate coefficient
  RCONST(565) = (2.50D-13*RO2)
  RCONST(566) = (8.06D-11)
  RCONST(567) = (SUN*3.98E-6)
  RCONST(568) = (3.00D-13)
  RCONST(569) = (2.20D-19)
  RCONST(570) = (4.45D-11)
  RCONST(571) = (SUN*2.68E-4)
  RCONST(572) = (1.40D-18*TEMP**2*EXP(194/TEMP))
  RCONST(573) = (1.4D-12*EXP(-1860/TEMP)*2.75)
  RCONST(574) = (2.00D-18*0.5)
  RCONST(575) = (2.00D-18*0.5)
  RCONST(576) = (6.2D-11*0.48)
  RCONST(577) = (6.2D-11*0.52)
  RCONST(578) = (SUN*8.92E-3*0.2*0.6)
  RCONST(579) = (SUN*8.92E-3*0.2*0.4)
  RCONST(580) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(581) = (5.00D-18*0.5)
  RCONST(582) = (5.00D-18*0.5)
  RCONST(583) = (7.10D-11*0.28)
  RCONST(584) = (7.10D-11*0.72)
  RCONST(585) = (SUN*8.92E-3*0.2*0.6)
  RCONST(586) = (SUN*8.92E-3*0.2*0.4)
  RCONST(587) = (1.00D-12)
  RCONST(588) = (1.40D-18)
  RCONST(589) = (3.01D-11)
  RCONST(590) = (1.00D-12)
  RCONST(591) = (8.00D-19)
  RCONST(592) = (6.90D-11)
  RCONST(593) = (3.98D-11)
  RCONST(594) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(595) = (4.28D-11)
  RCONST(596) = (SUN*5.79E-6+SUN*1.10E-6+SUN*1.10E-6)
  RCONST(597) = (3.96D-11)
  RCONST(598) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(599) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(600) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(601) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(602) = (7.5D-12*EXP(290/TEMP))
  RCONST(603) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(604) = (2.3E-12*1.74)
  RCONST(605) = (1.00D-11*RO2*0.70)
  RCONST(606) = (1.00D-11*RO2*0.30)
  RCONST(607) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(608) = (2.7D-12*EXP(360/TEMP))
! RCONST(609) = constant rate coefficient
  RCONST(610) = (9.20D-14*RO2)
  RCONST(611) = (1.2D-15)
  RCONST(612) = (1.0D-14)
  RCONST(613) = (1.0D-15)
  RCONST(614) = (7.0D-14)
  RCONST(615) = (6.0D-18)
  RCONST(616) = (1.0D-17)
  RCONST(617) = (9.33D-11)
  RCONST(618) = (SUN*5.79E-6+SUN*5.50E-5*2)
  RCONST(619) = (3.27D-11)
  RCONST(620) = (SUN*5.50E-5*2)
  RCONST(621) = (1.00D+06*0.5)
  RCONST(622) = (1.00D+06*0.5)
  RCONST(623) = (4.84D-11)
  RCONST(624) = (SUN*5.50E-5+SUN*1.25E-4)
  RCONST(625) = (8.11D-11)
  RCONST(626) = (SUN*5.50E-5*2)
  RCONST(627) = (4.28D-11)
  RCONST(628) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(629) = (SUN*3.98E-6)
  RCONST(630) = (SUN*5.79E-6)
  RCONST(631) = (1.00D+06)
  RCONST(632) = (5.24D-11)
  RCONST(633) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(634) = (SUN*3.98E-6)
  RCONST(635) = (1.43D-10)
  RCONST(636) = (SUN*5.79E-6)
  RCONST(637) = (SUN*2.99E-6)
  RCONST(638) = (1.00D+06)
  RCONST(639) = (1.40D-10)
  RCONST(640) = (SUN*2.99E-6)
  RCONST(641) = (3.13D-13)
  RCONST(642) = (1.53D-12)
  RCONST(643) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(644) = (2.7D-12*EXP(360/TEMP))
! RCONST(645) = constant rate coefficient
  RCONST(646) = (2.50D-13*RO2)
  RCONST(647) = (1.51D-10)
  RCONST(648) = (SUN*5.79E-6)
  RCONST(649) = (1.00D+06*0.52)
  RCONST(650) = (1.00D+06*0.48)
  RCONST(651) = (1.48D-10)
  RCONST(652) = (2.08D-12)
  RCONST(653) = (2.86D-13)
  RCONST(654) = (1.00D+06)
  RCONST(655) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(656) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(657) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(658) = (7.5D-12*EXP(290/TEMP))
  RCONST(659) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(660) = (2.3E-12*1.74)
  RCONST(661) = (1.00D-11*0.70*RO2)
  RCONST(662) = (1.00D-11*0.30*RO2)
  RCONST(663) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(664) = (2.7D-12*EXP(360/TEMP))
! RCONST(665) = constant rate coefficient
  RCONST(666) = (2.50D-13*RO2)
  RCONST(667) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(668) = (2.0D-18*0.5)
  RCONST(669) = (2.0D-18*0.5)
  RCONST(670) = (7.00D-11*0.5)
  RCONST(671) = (7.00D-11*0.5)
  RCONST(672) = (SUN*8.92E-3*0.2*0.6)
  RCONST(673) = (SUN*8.92E-3*0.2*0.4)
  RCONST(674) = (1.00D-12)
  RCONST(675) = (8.00D-19)
  RCONST(676) = (2.67D-11)
  RCONST(677) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(678) = (1.4D-12*EXP(-1860/TEMP)*4.25)
  RCONST(679) = (5.00D-18*0.50)
  RCONST(680) = (5.00D-18*0.50)
  RCONST(681) = (4.41D-11*0.385)
  RCONST(682) = (4.41D-11*0.23)
  RCONST(683) = (4.41D-11*0.385)
  RCONST(684) = (SUN*8.92E-3*0.2*0.3)
  RCONST(685) = (SUN*8.92E-3*0.2*0.3)
  RCONST(686) = (SUN*8.92E-3*0.2*0.4)
  RCONST(687) = (1.00D-12)
  RCONST(688) = (8.00D-19)
  RCONST(689) = (2.42D-11)
  RCONST(690) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(691) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(692) = (7.5D-12*EXP(290/TEMP))
  RCONST(693) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(694) = (2.3E-12*1.74)
  RCONST(695) = (1.00D-11*RO2)
  RCONST(696) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(697) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(698) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(699) = (7.5D-12*EXP(290/TEMP))
  RCONST(700) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(701) = (2.3E-12*1.74)
  RCONST(702) = (1.00D-11*0.70*RO2)
  RCONST(703) = (1.00D-11*0.30*RO2)
  RCONST(704) = (6.00D-11)
  RCONST(705) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(706) = (6.31D-11)
  RCONST(707) = (SUN*5.79E-6+SUN*1.80E-6+SUN*1.80E-6)
  RCONST(708) = (5.99D-11)
  RCONST(709) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(710) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(711) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(712) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(713) = (7.5D-12*EXP(290/TEMP))
  RCONST(714) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(715) = (2.3E-12*1.74)
  RCONST(716) = (1.00D-11*RO2*0.70)
  RCONST(717) = (1.00D-11*RO2*0.30)
  RCONST(718) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(719) = (2.7D-12*EXP(360/TEMP))
! RCONST(720) = constant rate coefficient
  RCONST(721) = (9.20D-14*RO2*0.7)
  RCONST(722) = (9.20D-14*RO2*0.3)
  RCONST(723) = (1.2D-15)
  RCONST(724) = (1.0D-14)
  RCONST(725) = (1.0D-15)
  RCONST(726) = (7.0D-14)
  RCONST(727) = (1.0D-17)
  RCONST(728) = (6.0D-18)
  RCONST(729) = (7.41D-11)
  RCONST(730) = (SUN*5.79E-6+SUN*5.50E-5+SUN*3.98E-6)
  RCONST(731) = (3.08D-11)
  RCONST(732) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(733) = (1.00D+06*0.5)
  RCONST(734) = (1.00D+06*0.5)
  RCONST(735) = (4.08D-11)
  RCONST(736) = (SUN*5.50E-5+SUN*2.68E-4)
  RCONST(737) = (6.20D-11)
  RCONST(738) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(739) = (7.36D-11)
  RCONST(740) = (SUN*5.79E-6+SUN*1.80E-6+SUN*1.80E-6)
  RCONST(741) = (1.00D+06)
  RCONST(742) = (5.22D-11)
  RCONST(743) = (SUN*2.68E-4+SUN*1.80E-6+SUN*1.80E-6)
  RCONST(744) = (9.76D-11)
  RCONST(745) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(746) = (GCARR(3.15E-14,0.0E+00,920.0))
  RCONST(747) = (3.13D-13)
  RCONST(748) = (1.53D-12)
  RCONST(749) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(750) = (2.7D-12*EXP(360/TEMP))
! RCONST(751) = constant rate coefficient
  RCONST(752) = (2.50D-13*RO2)
  RCONST(753) = (1.13D-10)
  RCONST(754) = (SUN*5.79E-6)
  RCONST(755) = (1.00D+06*0.85)
  RCONST(756) = (1.00D+06*0.15)
  RCONST(757) = (1.10D-10)
  RCONST(758) = (1.07D-10)
  RCONST(759) = (SUN*2.99E-6)
  RCONST(760) = (SUN*5.79E-6)
  RCONST(761) = (1.00D+06)
  RCONST(762) = (1.03D-10)
  RCONST(763) = (SUN*2.99E-6)
  RCONST(764) = (2.08D-12)
  RCONST(765) = (2.86D-13)
  RCONST(766) = (1.00D+06)
  RCONST(767) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(768) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(769) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(770) = (7.5D-12*EXP(290/TEMP))
  RCONST(771) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(772) = (2.3E-12*1.74)
  RCONST(773) = (1.00D-11*0.70*RO2)
  RCONST(774) = (1.00D-11*0.30*RO2)
  RCONST(775) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(776) = (2.7D-12*EXP(360/TEMP))
! RCONST(777) = constant rate coefficient
  RCONST(778) = (2.50D-13*RO2)
  RCONST(779) = (5.00D-18)
  RCONST(780) = (4.90D-11)
  RCONST(781) = (SUN*8.92E-3*0.0057)
  RCONST(782) = (5.87D-11)
  RCONST(783) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(784) = (6.18D-11)
  RCONST(785) = (SUN*5.79E-6+SUN*1.80E-6+SUN*1.80E-6)
  RCONST(786) = (5.85D-11)
  RCONST(787) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(788) = (7.26D-11)
  RCONST(789) = (SUN*5.79E-6+SUN*5.50E-5+SUN*3.98E-6)
  RCONST(790) = (2.92D-11)
  RCONST(791) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(792) = (1.00D+06*0.5)
  RCONST(793) = (1.00D+06*0.5)
  RCONST(794) = (3.58D-11)
  RCONST(795) = (SUN*5.50E-5+SUN*2.68E-4)
  RCONST(796) = (6.03D-11)
  RCONST(797) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(798) = (8.25D-11)
  RCONST(799) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(800) = (SUN*5.79E-6)
  RCONST(801) = (SUN*5.50E-5)
  RCONST(802) = (1.00D+06)
  RCONST(803) = (1.14D-10)
  RCONST(804) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(805) = (SUN*5.50E-5)
  RCONST(806) = (1.43D-10)
  RCONST(807) = (SUN*5.79E-6)
  RCONST(808) = (SUN*2.99E-6)
  RCONST(809) = (1.00D+06)
  RCONST(810) = (1.40D-10)
  RCONST(811) = (SUN*2.99E-6)
  RCONST(812) = (3.13D-13)
  RCONST(813) = (1.53D-12)
  RCONST(814) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(815) = (2.7D-12*EXP(360/TEMP))
! RCONST(816) = constant rate coefficient
  RCONST(817) = (2.50D-13*RO2)
  RCONST(818) = (1.51D-10)
  RCONST(819) = (SUN*5.79E-6)
  RCONST(820) = (1.00D+06*0.68)
  RCONST(821) = (1.00D+06*0.32)
  RCONST(822) = (1.48D-10)
  RCONST(823) = (2.08D-12)
  RCONST(824) = (2.86D-13)
  RCONST(825) = (1.00D+06)
  RCONST(826) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(827) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(828) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(829) = (7.5D-12*EXP(290/TEMP))
  RCONST(830) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(831) = (2.3E-12*1.74)
  RCONST(832) = (1.00D-11*0.30*RO2)
  RCONST(833) = (1.00D-11*0.70*RO2)
  RCONST(834) = (4.06D-11)
  RCONST(835) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(836) = (4.37D-11)
  RCONST(837) = (SUN*5.79E-6+SUN*1.10E-6+SUN*1.10E-6)
  RCONST(838) = (4.05D-11)
  RCONST(839) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(840) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(841) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(842) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(843) = (7.5D-12*EXP(290/TEMP))
  RCONST(844) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(845) = (2.3E-12*1.74)
  RCONST(846) = (1.00D-11*RO2*0.7)
  RCONST(847) = (1.00D-11*RO2*0.3)
  RCONST(848) = (2.91D-13*EXP(1300/TEMP)*0.520)
  RCONST(849) = (2.7D-12*EXP(360/TEMP))
! RCONST(850) = constant rate coefficient
  RCONST(851) = (8.80D-13*RO2*0.2)
  RCONST(852) = (8.80D-13*RO2*0.2)
  RCONST(853) = (8.80D-13*RO2*0.6)
  RCONST(854) = (1.31D-10)
  RCONST(855) = (SUN*5.79E-6+SUN*1.99E-5*2)
  RCONST(856) = (4.38D-11)
  RCONST(857) = (SUN*5.50E-5)
  RCONST(858) = (1.00D+06*0.5)
  RCONST(859) = (1.00D+06*0.5)
  RCONST(860) = (9.20D-11)
  RCONST(861) = (SUN*1.99E-5*2+SUN*3.98E-6)
  RCONST(862) = (8.23D-11)
  RCONST(863) = (SUN*1.99E-5*2)
  RCONST(864) = (7.52D-11)
  RCONST(865) = (SUN*5.79E-6)
  RCONST(866) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(867) = (1.00D+06)
  RCONST(868) = (4.90D-11)
  RCONST(869) = (SUN*1.25E-4+SUN*1.10E-6+SUN*1.10E-6)
  RCONST(870) = (7.75D-11)
  RCONST(871) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(872) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(873) = (2.7D-12*EXP(360/TEMP))
! RCONST(874) = constant rate coefficient
  RCONST(875) = (8.80D-13*0.20*RO2)
  RCONST(876) = (8.80D-13*0.20*RO2)
  RCONST(877) = (8.80D-13*0.60*RO2)
  RCONST(878) = (9.00D-14)
  RCONST(879) = (9.00D-13)
  RCONST(880) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(881) = (2.7D-12*EXP(360/TEMP))
! RCONST(882) = constant rate coefficient
  RCONST(883) = (2.50D-13*RO2)
  RCONST(884) = (1.07D-10)
  RCONST(885) = (SUN*2.99E-6)
  RCONST(886) = (SUN*5.79E-6)
  RCONST(887) = (1.00D+06)
  RCONST(888) = (1.04D-10)
  RCONST(889) = (SUN*2.99E-6)
  RCONST(890) = (2.08D-12)
  RCONST(891) = (2.86D-13)
  RCONST(892) = (1.00D+06)
  RCONST(893) = (1.16D-10)
  RCONST(894) = (SUN*5.79E-6)
  RCONST(895) = (1.00D+06*0.71)
  RCONST(896) = (1.00D+06*0.29)
  RCONST(897) = (1.13D-10)
  RCONST(898) = (7.5D-12*EXP(290/TEMP))
  RCONST(899) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(900) = (2.3E-12*1.74)
  RCONST(901) = (1.00D-11*RO2*0.7)
  RCONST(902) = (1.00D-11*RO2*0.3)
  RCONST(903) = (1.15D-10)
  RCONST(904) = (SUN*5.79E-6)
  RCONST(905) = (1.00D+06*0.68)
  RCONST(906) = (1.00D+06*0.32)
  RCONST(907) = (1.11D-10)
  RCONST(908) = (1.07D-10)
  RCONST(909) = (SUN*2.99E-6)
  RCONST(910) = (SUN*5.79E-6)
  RCONST(911) = (1.00D+06)
  RCONST(912) = (1.04D-10)
  RCONST(913) = (SUN*2.99E-6)
  RCONST(914) = (3.13D-13*1.0)
  RCONST(915) = (2.8D-12)
  RCONST(916) = (2.08D-12)
  RCONST(917) = (2.86D-13)
  RCONST(918) = (1.00D+06*1.0)
  RCONST(919) = (5.98D-11)
  RCONST(920) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(921) = (6.29D-11)
  RCONST(922) = (SUN*5.79E-6+SUN*1.10E-6+SUN*1.10E-6)
  RCONST(923) = (5.96D-11)
  RCONST(924) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(925) = (7.04D-11)
  RCONST(926) = (SUN*5.79E-6+SUN*5.50E-5+SUN*3.98E-6)
  RCONST(927) = (3.06D-11)
  RCONST(928) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(929) = (1.00D+06*0.50)
  RCONST(930) = (1.00D+06*0.50)
  RCONST(931) = (4.06D-11)
  RCONST(932) = (SUN*5.50E-5+SUN*2.68E-4)
  RCONST(933) = (6.17D-11)
  RCONST(934) = (SUN*5.50E-5+SUN*3.98E-6)
  RCONST(935) = (9.42D-11)
  RCONST(936) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(937) = (SUN*5.79E-6)
  RCONST(938) = (1.00D+06)
  RCONST(939) = (6.45D-11)
  RCONST(940) = (SUN*1.25E-4+SUN*1.80E-6+SUN*1.80E-6)
  RCONST(941) = (1.18D-10)
  RCONST(942) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(943) = (1.10D-12)
  RCONST(944) = (4.66D-12)
  RCONST(945) = (SUN*5.79E-6)
  RCONST(946) = (1.06D-12)
  RCONST(947) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3)*0.67)
  RCONST(948) = (4.65D-11)
  RCONST(949) = (SUN*5.79E-6)
  RCONST(950) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(951) = (2.7D-12*EXP(360/TEMP))
! RCONST(952) = constant rate coefficient
  RCONST(953) = (8.80D-13*RO2)
  RCONST(954) = (1.00D+06*0.5)
  RCONST(955) = (1.00D+06*0.5)
  RCONST(956) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(957) = (2.7D-12*EXP(360/TEMP))
! RCONST(958) = constant rate coefficient
  RCONST(959) = (8.80D-13*RO2*0.2)
  RCONST(960) = (8.80D-13*RO2*0.6)
  RCONST(961) = (8.80D-13*RO2*0.2)
  RCONST(962) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(963) = (2.7D-12*EXP(360/TEMP))
! RCONST(964) = constant rate coefficient
  RCONST(965) = (2.00D-12*0.6*RO2)
  RCONST(966) = (2.00D-12*0.2*RO2)
  RCONST(967) = (2.00D-12*0.2*RO2)
  RCONST(968) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(969) = (2.7D-12*EXP(360/TEMP))
! RCONST(970) = constant rate coefficient
  RCONST(971) = (1.40D-12*0.20*RO2)
  RCONST(972) = (1.40D-12*0.20*RO2)
  RCONST(973) = (1.40D-12*0.60*RO2)
  RCONST(974) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(975) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(976) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(977) = (7.5D-12*EXP(290/TEMP))
  RCONST(978) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(979) = (2.3E-12*1.74)
  RCONST(980) = (1.00D-11*0.30*RO2)
  RCONST(981) = (1.00D-11*0.70*RO2)
  RCONST(982) = (1.00D+06*0.11)
  RCONST(983) = (1.00D+06*0.89)
  RCONST(984) = (1.00D+06*0.27)
  RCONST(985) = (1.00D+06*0.255)
  RCONST(986) = (1.00D+06*0.255)
  RCONST(987) = (1.00D+06*0.22)
  RCONST(988) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(989) = (2.7D-12*EXP(360/TEMP))
! RCONST(990) = constant rate coefficient
  RCONST(991) = (9.20D-14*0.3*RO2)
  RCONST(992) = (9.20D-14*0.7*RO2)
  RCONST(993) = (1.00D-12)
  RCONST(994) = (8.00D-17)
  RCONST(995) = (2.66D-11)
  RCONST(996) = (7.00D11*EXP(-3160/TEMP)+5.00D-12)
  RCONST(997) = (5.00D-12*3.2*(1-EXP(-550/TEMP)))
  RCONST(998) = (5.00D-12*3.2*EXP(-550/TEMP))
  RCONST(999) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1000) = (2.7D-12*EXP(360/TEMP))
! RCONST(1001) = constant rate coefficient
  RCONST(1002) = (9.20D-14*RO2)
  RCONST(1003) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1004) = (2.7D-12*EXP(360/TEMP))
  RCONST(1005) = (2.3E-12*1.74)
  RCONST(1006) = (8.80D-13*RO2*0.2)
  RCONST(1007) = (8.80D-13*RO2*0.6)
  RCONST(1008) = (8.80D-13*RO2*0.2)
  RCONST(1009) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1010) = (2.7D-12*EXP(360/TEMP))
! RCONST(1011) = constant rate coefficient
  RCONST(1012) = (9.20D-14*RO2*0.70)
  RCONST(1013) = (9.20D-14*RO2*0.30)
  RCONST(1014) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1015) = (2.7D-12*EXP(360/TEMP))
! RCONST(1016) = constant rate coefficient
  RCONST(1017) = (9.20D-14*RO2)
  RCONST(1018) = (1.00D+06*0.5)
  RCONST(1019) = (1.00D+06*0.5)
  RCONST(1020) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1021) = (2.7D-12*EXP(360/TEMP))
! RCONST(1022) = constant rate coefficient
  RCONST(1023) = (9.20D-14*RO2*0.7)
  RCONST(1024) = (9.20D-14*RO2*0.3)
  RCONST(1025) = (2.22D-11)
  RCONST(1026) = (SUN*5.50E-5)
  RCONST(1027) = (2.53D-11)
  RCONST(1028) = (SUN*5.79E-6+SUN*5.50E-5)
  RCONST(1029) = (2.21D-11)
  RCONST(1030) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1031) = (6.97D-11)
  RCONST(1032) = (SUN*5.50E-5)
  RCONST(1033) = (SUN*5.79E-6)
  RCONST(1034) = (SUN*3.98E-6)
  RCONST(1035) = (1.00D+06*0.5)
  RCONST(1036) = (1.00D+06*0.5)
  RCONST(1037) = (SUN*1.25E-4)
  RCONST(1038) = (1.23D-11)
  RCONST(1039) = (1.36D-10)
  RCONST(1040) = (SUN*1.99E-5*2)
  RCONST(1041) = (6.63D-11)
  RCONST(1042) = (SUN*5.50E-5)
  RCONST(1043) = (SUN*3.98E-6)
  RCONST(1044) = (5.08D-11)
  RCONST(1045) = (SUN*1.80E-6)
  RCONST(1046) = (2.08D-12)
  RCONST(1047) = (2.86D-13)
  RCONST(1048) = (4.65D-11)
  RCONST(1049) = (SUN*5.79E-6)
  RCONST(1050) = (1.00D-12)
  RCONST(1051) = (2.35D-10)
  RCONST(1052) = (5.03D-12)
  RCONST(1053) = (6.83D-12)
  RCONST(1054) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1055) = (2.7D-12*EXP(360/TEMP))
! RCONST(1056) = constant rate coefficient
  RCONST(1057) = (8.80D-13*RO2)
  RCONST(1058) = (5.44D-11)
  RCONST(1059) = (SUN*1.80E-6+SUN*1.80E-6)
  RCONST(1060) = (SUN*1.25E-4)
  RCONST(1061) = (8.0D-13)
  RCONST(1062) = (1.10D-12)
  RCONST(1063) = (4.66D-12)
  RCONST(1064) = (SUN*5.79E-6)
  RCONST(1065) = (1.06D-12)
  RCONST(1066) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3)*0.67)
  RCONST(1067) = (1.90D-12*EXP(190/TEMP))
  RCONST(1068) = (SUN*5.79E-6)
  RCONST(1069) = (3.90D-13)
  RCONST(1070) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1071) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1072) = (7.5D-12*EXP(290/TEMP))
  RCONST(1073) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1074) = (2.3E-12*1.74)
  RCONST(1075) = (1.00D-11*0.70*RO2)
  RCONST(1076) = (1.00D-11*0.30*RO2)
  RCONST(1077) = (1.00D+06*0.89)
  RCONST(1078) = (1.00D+06*0.11)
  RCONST(1079) = (1.00D+06*0.11)
  RCONST(1080) = (1.00D+06*0.89)
  RCONST(1081) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1082) = (2.7D-12*EXP(360/TEMP))
! RCONST(1083) = constant rate coefficient
  RCONST(1084) = (9.20D-14*0.70*RO2)
  RCONST(1085) = (9.20D-14*0.30*RO2)
  RCONST(1086) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1087) = (2.7D-12*EXP(360/TEMP))
! RCONST(1088) = constant rate coefficient
  RCONST(1089) = (9.20D-14*RO2)
  RCONST(1090) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1091) = (2.7D-12*EXP(360/TEMP))
! RCONST(1092) = constant rate coefficient
  RCONST(1093) = (8.80D-13*RO2*0.2)
  RCONST(1094) = (8.80D-13*RO2*0.6)
  RCONST(1095) = (8.80D-13*RO2*0.2)
  RCONST(1096) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1097) = (2.7D-12*EXP(360/TEMP))
! RCONST(1098) = constant rate coefficient
  RCONST(1099) = (9.20D-14*RO2*0.70)
  RCONST(1100) = (9.20D-14*RO2*0.30)
  RCONST(1101) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(1102) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1103) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(1104) = (7.5D-12*EXP(290/TEMP))
  RCONST(1105) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1106) = (2.3E-12*1.74)
  RCONST(1107) = (1.00D-11*0.70*RO2)
  RCONST(1108) = (1.00D-11*0.30*RO2)
  RCONST(1109) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1110) = (2.7D-12*EXP(360/TEMP))
! RCONST(1111) = constant rate coefficient
  RCONST(1112) = (9.20D-14*RO2*0.30)
  RCONST(1113) = (9.20D-14*RO2*0.70)
  RCONST(1114) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1115) = (2.7D-12*EXP(360/TEMP))
! RCONST(1116) = constant rate coefficient
  RCONST(1117) = (9.20D-14*RO2)
  RCONST(1118) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(1119) = (2.7D-12*EXP(360/TEMP))
! RCONST(1120) = constant rate coefficient
  RCONST(1121) = (2.00D-12*RO2)
  RCONST(1122) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1123) = (2.7D-12*EXP(360/TEMP))
! RCONST(1124) = constant rate coefficient
  RCONST(1125) = (9.20D-14*RO2*0.70)
  RCONST(1126) = (9.20D-14*RO2*0.30)
  RCONST(1127) = (1.00D+06*0.35)
  RCONST(1128) = (1.00D+06*0.65)
  RCONST(1129) = (2.40D-17*0.27)
  RCONST(1130) = (2.40D-17*0.73)
  RCONST(1131) = (4.72D-11)
  RCONST(1132) = (SUN*5.64E-4*2)
  RCONST(1133) = (4.37D-11)
  RCONST(1134) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1135) = (1.00D+06*0.83)
  RCONST(1136) = (1.00D+06*0.17)
  RCONST(1137) = (3.59D-12)
  RCONST(1138) = (SUN*5.64E-4*2)
  RCONST(1139) = (5.43D-11)
  RCONST(1140) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1141) = (2.27D-11)
  RCONST(1142) = (SUN*5.50E-5)
  RCONST(1143) = (2.57D-11)
  RCONST(1144) = (SUN*5.79E-6+SUN*5.50E-5)
  RCONST(1145) = (2.25D-11)
  RCONST(1146) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1147) = (1.35D-10)
  RCONST(1148) = (SUN*5.79E-6)
  RCONST(1149) = (SUN*5.50E-5*2)
  RCONST(1150) = (1.00D+06)
  RCONST(1151) = (4.18D-11)
  RCONST(1152) = (SUN*5.50E-5*2)
  RCONST(1153) = (2.91D-13*EXP(1300/TEMP)*0.387)
  RCONST(1154) = (2.7D-12*EXP(360/TEMP))
! RCONST(1155) = constant rate coefficient
  RCONST(1156) = (2.00D-12*0.2*RO2)
  RCONST(1157) = (2.00D-12*0.6*RO2)
  RCONST(1158) = (2.00D-12*0.2*RO2)
  RCONST(1159) = (SUN*1.99E-5)
  RCONST(1160) = (1.4D-12*EXP(-1860/TEMP)*4.0)
  RCONST(1161) = (2.45D-11)
  RCONST(1162) = (2.08D-12)
  RCONST(1163) = (2.86D-13)
  RCONST(1164) = (4.65D-11)
  RCONST(1165) = (SUN*5.79E-6)
  RCONST(1166) = (1.00D-12)
  RCONST(1167) = (2.35D-11)
  RCONST(1168) = (5.03D-12)
  RCONST(1169) = (6.83D-12)
  RCONST(1170) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1171) = (2.7D-12*EXP(360/TEMP))
! RCONST(1172) = constant rate coefficient
  RCONST(1173) = (8.80D-13*RO2)
  RCONST(1174) = (4.66D-11)
  RCONST(1175) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(1176) = (1.10D-12)
  RCONST(1177) = (4.66D-12)
  RCONST(1178) = (SUN*5.79E-6)
  RCONST(1179) = (1.06D-12)
  RCONST(1180) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3)*0.67)
  RCONST(1181) = (1.90D-12*EXP(190/TEMP))
  RCONST(1182) = (SUN*5.79E-6)
  RCONST(1183) = (3.90D-13)
  RCONST(1184) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1185) = (2.7D-12*EXP(360/TEMP))
! RCONST(1186) = constant rate coefficient
  RCONST(1187) = (8.80D-13*0.20*RO2)
  RCONST(1188) = (8.80D-13*0.20*RO2)
  RCONST(1189) = (8.80D-13*0.60*RO2)
  RCONST(1190) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(1191) = (2.7D-12*EXP(360/TEMP))
! RCONST(1192) = constant rate coefficient
  RCONST(1193) = (8.80D-12*RO2)
  RCONST(1194) = (2.08D-12)
  RCONST(1195) = (2.86D-13)
  RCONST(1196) = (4.65D-11)
  RCONST(1197) = (SUN*5.79E-6)
  RCONST(1198) = (1.00D-12)
  RCONST(1199) = (2.35D-11)
  RCONST(1200) = (5.03D-12)
  RCONST(1201) = (6.83D-12)
  RCONST(1202) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1203) = (2.7D-12*EXP(360/TEMP))
! RCONST(1204) = constant rate coefficient
  RCONST(1205) = (8.80D-13*RO2)
  RCONST(1206) = (4.38D-11)
  RCONST(1207) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(1208) = (3.70D-11)
  RCONST(1209) = (SUN*1.10E-6+SUN*1.10E-6)
  RCONST(1210) = (4.00D-11)
  RCONST(1211) = (SUN*5.64E-4*2)
  RCONST(1212) = (1.00D+06*0.40)
  RCONST(1213) = (1.00D+06*0.60)
  RCONST(1214) = (3.70D-11)
  RCONST(1215) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1216) = (2.31D-11)
  RCONST(1217) = (SUN*5.50E-5)
  RCONST(1218) = (2.62D-11)
  RCONST(1219) = (SUN*5.79E-6+SUN*5.50E-5)
  RCONST(1220) = (2.29D-11)
  RCONST(1221) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1222) = (1.44D-10)
  RCONST(1223) = (SUN*5.79E-6)
  RCONST(1224) = (SUN*1.99E-5*2)
  RCONST(1225) = (1.00D+06)
  RCONST(1226) = (5.77D-11)
  RCONST(1227) = (SUN*1.99E-5*2)
  RCONST(1228) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1229) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1230) = (7.5D-12*EXP(290/TEMP))
  RCONST(1231) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1232) = (2.3E-12*1.74)
  RCONST(1233) = (1.00D-11*RO2)
  RCONST(1234) = (1.22D-10)
  RCONST(1235) = (1.90D-12*EXP(190/TEMP))
  RCONST(1236) = (SUN*1.99E-5*2)
  RCONST(1237) = (SUN*5.79E-6)
  RCONST(1238) = (1.00D+06)
  RCONST(1239) = (3.67D-11)
  RCONST(1240) = (SUN*1.25E-4+SUN*1.99E-5)
  RCONST(1241) = (8.13D-11)
  RCONST(1242) = (SUN*1.99E-5*2)
  RCONST(1243) = (2.08D-12)
  RCONST(1244) = (2.86D-13)
  RCONST(1245) = (3.60D-12)
  RCONST(1246) = (SUN*5.79E-6)
  RCONST(1247) = (2.60D-12)
  RCONST(1248) = (3.47D-12)
  RCONST(1249) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1250) = (2.7D-12*EXP(360/TEMP))
! RCONST(1251) = constant rate coefficient
  RCONST(1252) = (8.80D-13*RO2)
  RCONST(1253) = (3.00D-13)
  RCONST(1254) = (4.6D-12)
  RCONST(1255) = (4.66D-12)
  RCONST(1256) = (SUN*5.79E-6)
  RCONST(1257) = (1.10D-12)
  RCONST(1258) = (1.06D-12)
  RCONST(1259) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3)*0.67)
  RCONST(1260) = (1.00D-12)
  RCONST(1261) = (2.3D-11)
  RCONST(1262) = (2.08D-12)
  RCONST(1263) = (2.86D-13)
  RCONST(1264) = (5.03D-12)
  RCONST(1265) = (6.83D-12)
  RCONST(1266) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1267) = (2.7D-12*EXP(360/TEMP))
! RCONST(1268) = constant rate coefficient
  RCONST(1269) = (8.80D-13*RO2)
  RCONST(1270) = (6.18D-12)
  RCONST(1271) = (SUN*5.79E-6)
  RCONST(1272) = (1.00D+06*0.5)
  RCONST(1273) = (1.00D+06*0.5)
  RCONST(1274) = (1.20D-15)
  RCONST(1275) = (1.00D-14)
  RCONST(1276) = (1.00D-15)
  RCONST(1277) = (7.00D-14)
  RCONST(1278) = (6.00D-18)
  RCONST(1279) = (1.00D-17)
  RCONST(1280) = (3.68D-11)
  RCONST(1281) = (SUN*5.79E-6)
  RCONST(1282) = (1.00D+06)
  RCONST(1283) = (1.78D-11)
  RCONST(1284) = (2.78D-11)
  RCONST(1285) = (SUN*5.79E-6+SUN*2.68E-4)
  RCONST(1286) = (1.90D-12*EXP(190/TEMP))
  RCONST(1287) = (5.99D-12)
  RCONST(1288) = (1.00D+06)
  RCONST(1289) = (SUN*2.68E-4)
  RCONST(1290) = (2.69D-12)
  RCONST(1291) = (SUN*1.25E-4)
  RCONST(1292) = (SUN*2.68E-4)
  RCONST(1293) = (1.4D-12*EXP(-1860/TEMP)*4.0)
  RCONST(1294) = (1.23D-11)
  RCONST(1295) = (1.21D-10)
  RCONST(1296) = (1.90D-12*EXP(190/TEMP))
  RCONST(1297) = (SUN*5.79E-6+SUN*1.99E-5+SUN*3.98E-6)
  RCONST(1298) = (1.00D+06)
  RCONST(1299) = (7.48D-11)
  RCONST(1300) = (SUN*1.99E-5+SUN*2.68E-4)
  RCONST(1301) = (6.26D-11)
  RCONST(1302) = (SUN*1.99E-5+SUN*3.98E-6)
  RCONST(1303) = (5.44D-11)
  RCONST(1304) = (SUN*5.64E-4*2)
  RCONST(1305) = (1.00D+06*0.35)
  RCONST(1306) = (1.00D+06*0.65)
  RCONST(1307) = (5.08D-11)
  RCONST(1308) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1309) = (8.8D-12*EXP(-1320/TEMP)+1.7D-14*EXP(423/TEMP))
  RCONST(1310) = (SUN*6.09E-7)
  RCONST(1311) = (1.20D-15)
  RCONST(1312) = (1.00D-14)
  RCONST(1313) = (1.00D-15)
  RCONST(1314) = (7.00D-14)
  RCONST(1315) = (6.00D-18)
  RCONST(1316) = (3.12D-11)
  RCONST(1317) = (SUN*3.98E-6+SUN*5.50E-5)
  RCONST(1318) = (SUN*5.79E-6)
  RCONST(1319) = (1.00D+06)
  RCONST(1320) = (2.77D-11)
  RCONST(1321) = (SUN*3.98E-6+SUN*5.50E-5)
  RCONST(1322) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1323) = (2.7D-12*EXP(360/TEMP))
! RCONST(1324) = constant rate coefficient
  RCONST(1325) = (9.20D-14*RO2)
  RCONST(1326) = (1.85D-11)
  RCONST(1327) = (SUN*3.98E-6+SUN*1.25E-4)
  RCONST(1328) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1329) = (2.7D-12*EXP(360/TEMP))
! RCONST(1330) = constant rate coefficient
  RCONST(1331) = (9.20D-14*RO2*0.70)
  RCONST(1332) = (9.20D-14*RO2*0.30)
  RCONST(1333) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(1334) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(1335) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1336) = (7.5D-12*EXP(290/TEMP))
  RCONST(1337) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1338) = (2.3E-12*1.74)
  RCONST(1339) = (1.00D-11*0.7*RO2)
  RCONST(1340) = (1.00D-11*0.3*RO2)
  RCONST(1341) = (4.95D-12)
  RCONST(1342) = (SUN*5.79E-6)
  RCONST(1343) = (1.00D+06)
  RCONST(1344) = (4.42D-11)
  RCONST(1345) = (1.90D-12*EXP(190/TEMP))
  RCONST(1346) = (SUN*5.79E-6)
  RCONST(1347) = (1.00D+06)
  RCONST(1348) = (4.65D-13)
  RCONST(1349) = (SUN*3.98E-6)
  RCONST(1350) = (1.91D-11)
  RCONST(1351) = (7.12D-12)
  RCONST(1352) = (SUN*5.79E-6)
  RCONST(1353) = (1.00D+06)
  RCONST(1354) = (3.76D-12)
  RCONST(1355) = (4.44D-12)
  RCONST(1356) = (SUN*5.79E-6)
  RCONST(1357) = (1.00D+06)
  RCONST(1358) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(1359) = (2.7D-12*EXP(360/TEMP))
! RCONST(1360) = constant rate coefficient
  RCONST(1361) = (2.00D-12*RO2)
  RCONST(1362) = (1.20D-15)
  RCONST(1363) = (1.00D-14)
  RCONST(1364) = (1.00D-15)
  RCONST(1365) = (7.00D-14)
  RCONST(1366) = (6.00D-18)
  RCONST(1367) = (1.00D-17)
  RCONST(1368) = (2.53D-11)
  RCONST(1369) = (SUN*5.79E-6)
  RCONST(1370) = (1.00D+06)
  RCONST(1371) = (2.19D-11)
  RCONST(1372) = (7.83D-15)
  RCONST(1373) = (5.10D-14)
  RCONST(1374) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1375) = (2.7D-12*EXP(360/TEMP))
! RCONST(1376) = constant rate coefficient
  RCONST(1377) = (2.50D-13*RO2)
  RCONST(1378) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1379) = (2.7D-12*EXP(360/TEMP))
! RCONST(1380) = constant rate coefficient
  RCONST(1381) = (8.80D-13*RO2)
  RCONST(1382) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1383) = (2.7D-12*EXP(360/TEMP))
! RCONST(1384) = constant rate coefficient
  RCONST(1385) = (8.80D-13*RO2*0.2)
  RCONST(1386) = (8.80D-13*RO2*0.6)
  RCONST(1387) = (8.80D-13*RO2*0.2)
  RCONST(1388) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1389) = (2.7D-12*EXP(360/TEMP))
! RCONST(1390) = constant rate coefficient
  RCONST(1391) = (8.00D-13*RO2)
  RCONST(1392) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1393) = (2.7D-12*EXP(360/TEMP))
! RCONST(1394) = constant rate coefficient
  RCONST(1395) = (8.00D-13*RO2)
  RCONST(1396) = (2.05D-10)
  RCONST(1397) = (SUN*5.79E-6)
  RCONST(1398) = (1.00D+06*0.35)
  RCONST(1399) = (1.00D+06*0.65)
  RCONST(1400) = (3.90D-11)
  RCONST(1401) = (SUN*5.64E-4*2)
  RCONST(1402) = (3.54D-11)
  RCONST(1403) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1404) = (3.54D-11)
  RCONST(1405) = (SUN*1.80E-6)
  RCONST(1406) = (1.2D-15)
  RCONST(1407) = (1.0D-14)
  RCONST(1408) = (1.0D-15)
  RCONST(1409) = (7.0D-14)
  RCONST(1410) = (6.0D-18)
  RCONST(1411) = (4.46D-11)
  RCONST(1412) = (SUN*5.50E-5)
  RCONST(1413) = (SUN*5.79E-6+SUN*3.98E-6)
  RCONST(1414) = (1.00D+06)
  RCONST(1415) = (4.13D-11)
  RCONST(1416) = (SUN*3.98E-6+SUN*5.50E-5)
  RCONST(1417) = (5.80D-12)
  RCONST(1418) = (SUN*5.79E-6)
  RCONST(1419) = (1.00D+06)
  RCONST(1420) = (6.01D-11)
  RCONST(1421) = (1.90D-12*EXP(190/TEMP))
  RCONST(1422) = (SUN*5.79E-6)
  RCONST(1423) = (1.00D+06)
  RCONST(1424) = (1.09D-13)
  RCONST(1425) = (SUN*3.98E-6)
  RCONST(1426) = (2.54D-11)
  RCONST(1427) = (3.09D-11)
  RCONST(1428) = (SUN*5.79E-6)
  RCONST(1429) = (1.00D+06)
  RCONST(1430) = (2.74D-11)
  RCONST(1431) = (2.40D-17*0.73)
  RCONST(1432) = (2.40D-17*0.27)
  RCONST(1433) = (4.72D-11)
  RCONST(1434) = (SUN*5.64E-4*2)
  RCONST(1435) = (1.00D+06*0.35)
  RCONST(1436) = (1.00D+06*0.65)
  RCONST(1437) = (4.37D-11)
  RCONST(1438) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1439) = (4.58D-11)
  RCONST(1440) = (SUN*5.79E-6)
  RCONST(1441) = (SUN*5.50E-5*2)
  RCONST(1442) = (1.00D+06)
  RCONST(1443) = (5.70D-11)
  RCONST(1444) = (SUN*5.50E-5*2)
  RCONST(1445) = (5.16D-12)
  RCONST(1446) = (SUN*5.79E-6)
  RCONST(1447) = (1.00D+06)
  RCONST(1448) = (2.00D-11)
  RCONST(1449) = (SUN*3.98E-6)
  RCONST(1450) = (SUN*5.79E-6)
  RCONST(1451) = (1.00D+06)
  RCONST(1452) = (2.78D-11)
  RCONST(1453) = (SUN*5.79E-6)
  RCONST(1454) = (1.00D+06)
  RCONST(1455) = (2.44D-11)
  RCONST(1456) = (1.50D-12)
  RCONST(1457) = (3.69D-12)
  RCONST(1458) = (SUN*5.79E-6+SUN*2.68E-4)
  RCONST(1459) = (1.00D+06)
  RCONST(1460) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1461) = (2.7D-12*EXP(360/TEMP))
! RCONST(1462) = constant rate coefficient
  RCONST(1463) = (9.20D-14*RO2)
  RCONST(1464) = (1.4D-12)
  RCONST(1465) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1466) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1467) = (7.5D-12*EXP(290/TEMP))
  RCONST(1468) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1469) = (2.3E-12*1.74)
  RCONST(1470) = (1.00D-11*RO2)
  RCONST(1471) = (2.91D-11)
  RCONST(1472) = (1.90D-12*EXP(190/TEMP))
  RCONST(1473) = (SUN*5.79E-6)
  RCONST(1474) = (SUN*1.99E-5)
  RCONST(1475) = (1.00D+06)
  RCONST(1476) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1477) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1478) = (7.5D-12*EXP(290/TEMP))
  RCONST(1479) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1480) = (2.3E-12*1.74)
  RCONST(1481) = (1.00D-11*RO2)
  RCONST(1482) = (7.83D-15)
  RCONST(1483) = (5.10D-14)
  RCONST(1484) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1485) = (2.7D-12*EXP(360/TEMP))
! RCONST(1486) = constant rate coefficient
  RCONST(1487) = (2.50D-13*RO2)
  RCONST(1488) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1489) = (2.7D-12*EXP(360/TEMP))
! RCONST(1490) = constant rate coefficient
  RCONST(1491) = (9.20D-14*RO2)
  RCONST(1492) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1493) = (2.7D-12*EXP(360/TEMP))
! RCONST(1494) = constant rate coefficient
  RCONST(1495) = (9.20D-14*RO2*0.7)
  RCONST(1496) = (9.20D-14*RO2*0.3)
  RCONST(1497) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1498) = (2.7D-12*EXP(360/TEMP))
! RCONST(1499) = constant rate coefficient
  RCONST(1500) = (8.00D-13*RO2)
  RCONST(1501) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1502) = (2.7D-12*EXP(360/TEMP))
! RCONST(1503) = constant rate coefficient
  RCONST(1504) = (8.00D-13*RO2)
  RCONST(1505) = (2.05D-10)
  RCONST(1506) = (SUN*5.79E-6)
  RCONST(1507) = (3.6D-12)
  RCONST(1508) = (9.75D-11)
  RCONST(1509) = (1.90D-12*EXP(190/TEMP))
  RCONST(1510) = (SUN*3.98E-6*2)
  RCONST(1511) = (SUN*5.79E-6)
  RCONST(1512) = (1.00D+06*0.7)
  RCONST(1513) = (1.00D+06*0.3)
  RCONST(1514) = (4.00D-11)
  RCONST(1515) = (SUN*3.98E-6*2)
  RCONST(1516) = (9.04D-12)
  RCONST(1517) = (SUN*3.98E-6)
  RCONST(1518) = (SUN*2.68E-4)
  RCONST(1519) = (7.83D-11)
  RCONST(1520) = (SUN*5.79E-6)
  RCONST(1521) = (SUN*3.98E-6)
  RCONST(1522) = (SUN*5.50E-5)
  RCONST(1523) = (1.00D+06*0.5)
  RCONST(1524) = (1.00D+06*0.5)
  RCONST(1525) = (7.83D-15)
  RCONST(1526) = (5.10D-14)
  RCONST(1527) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1528) = (2.7D-12*EXP(360/TEMP))
! RCONST(1529) = constant rate coefficient
  RCONST(1530) = (2.50D-13*RO2)
  RCONST(1531) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1532) = (2.7D-12*EXP(360/TEMP))
! RCONST(1533) = constant rate coefficient
  RCONST(1534) = (9.20D-14*RO2)
  RCONST(1535) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1536) = (2.7D-12*EXP(360/TEMP))
! RCONST(1537) = constant rate coefficient
  RCONST(1538) = (9.20D-14*RO2*0.70)
  RCONST(1539) = (9.20D-14*RO2*0.30)
  RCONST(1540) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1541) = (2.7D-12*EXP(360/TEMP))
! RCONST(1542) = constant rate coefficient
  RCONST(1543) = (8.00D-13*RO2)
  RCONST(1544) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1545) = (2.7D-12*EXP(360/TEMP))
! RCONST(1546) = constant rate coefficient
  RCONST(1547) = (8.00D-13*RO2)
  RCONST(1548) = (2.05D-10)
  RCONST(1549) = (SUN*5.79E-6)
  RCONST(1550) = (7.33D-11)
  RCONST(1551) = (SUN*5.79E-6)
  RCONST(1552) = (6.97D-11)
  RCONST(1553) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1554) = (2.45D-11)
  RCONST(1555) = (SUN*1.25E-4*2)
  RCONST(1556) = (2.25D-15)
  RCONST(1557) = (3.00D-14)
  RCONST(1558) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1559) = (2.7D-12*EXP(360/TEMP))
! RCONST(1560) = constant rate coefficient
  RCONST(1561) = (2.50D-13*RO2)
  RCONST(1562) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1563) = (2.7D-12*EXP(360/TEMP))
! RCONST(1564) = constant rate coefficient
  RCONST(1565) = (8.00D-13*RO2)
  RCONST(1566) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1567) = (2.7D-12*EXP(360/TEMP))
! RCONST(1568) = constant rate coefficient
  RCONST(1569) = (8.00D-13*RO2)
  RCONST(1570) = (1.90D-12*EXP(190/TEMP))
  RCONST(1571) = (SUN*5.79E-6)
  RCONST(1572) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1573) = (2.7D-12*EXP(360/TEMP))
! RCONST(1574) = constant rate coefficient
  RCONST(1575) = (8.80D-13*RO2)
  RCONST(1576) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1577) = (2.7D-12*EXP(360/TEMP))
! RCONST(1578) = constant rate coefficient
  RCONST(1579) = (8.80D-13*RO2*0.2)
  RCONST(1580) = (8.80D-13*RO2*0.6)
  RCONST(1581) = (8.80D-13*RO2*0.2)
  RCONST(1582) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1583) = (2.7D-12*EXP(360/TEMP))
! RCONST(1584) = constant rate coefficient
  RCONST(1585) = (8.80D-13*RO2)
  RCONST(1586) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1587) = (2.7D-12*EXP(360/TEMP))
! RCONST(1588) = constant rate coefficient
  RCONST(1589) = (8.80D-13*RO2*0.2)
  RCONST(1590) = (8.80D-13*RO2*0.6)
  RCONST(1591) = (8.80D-13*RO2*0.2)
  RCONST(1592) = (7.83D-15)
  RCONST(1593) = (5.10D-14)
  RCONST(1594) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1595) = (2.7D-12*EXP(360/TEMP))
! RCONST(1596) = constant rate coefficient
  RCONST(1597) = (2.50D-13*RO2)
  RCONST(1598) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1599) = (2.7D-12*EXP(360/TEMP))
! RCONST(1600) = constant rate coefficient
  RCONST(1601) = (8.00D-13*RO2)
  RCONST(1602) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1603) = (2.7D-12*EXP(360/TEMP))
! RCONST(1604) = constant rate coefficient
  RCONST(1605) = (8.00D-13*RO2)
  RCONST(1606) = (2.05D-10)
  RCONST(1607) = (SUN*5.79E-6)
  RCONST(1608) = (1.4D-12*EXP(-1860/TEMP)*8.0)
  RCONST(1609) = (3.44D-11)
  RCONST(1610) = (SUN*1.25E-4)
  RCONST(1611) = (1.16D-12)
  RCONST(1612) = (2.19D-11)
  RCONST(1613) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1614) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(1615) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(1616) = (7.5D-12*EXP(290/TEMP))
  RCONST(1617) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1618) = (2.3E-12*1.74)
  RCONST(1619) = (1.00D-11*0.7*RO2)
  RCONST(1620) = (1.00D-11*0.3*RO2)
  RCONST(1621) = (1.4D-12*EXP(-1860/TEMP)*5.5)
  RCONST(1622) = (1.72D-11)
  RCONST(1623) = (SUN*1.25E-4+SUN*2.68E-4*2)
  RCONST(1624) = (1.36D-13*EXP(1250/TEMP)*0.15)
  RCONST(1625) = (1.36D-13*EXP(1250/TEMP)*0.85)
  RCONST(1626) = (2.7D-12*EXP(360/TEMP))
! RCONST(1627) = constant rate coefficient
  RCONST(1628) = (2*(3.5D-13*8.0D-12)**0.5*RO2*0.2)
  RCONST(1629) = (2*(3.5D-13*8.0D-12)**0.5*RO2*0.6)
  RCONST(1630) = (2*(3.5D-13*8.0D-12)**0.5*RO2*0.2)
  RCONST(1631) = (9.93D-11)
  RCONST(1632) = (SUN*5.79E-6)
  RCONST(1633) = (1.00D+06)
  RCONST(1634) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(1635) = (2.7D-12*EXP(360/TEMP))
! RCONST(1636) = constant rate coefficient
  RCONST(1637) = (8.80D-13*0.2*RO2)
  RCONST(1638) = (8.80D-13*0.2*RO2)
  RCONST(1639) = (8.80D-13*0.6*RO2)
  RCONST(1640) = (1.37D-11)
  RCONST(1641) = (SUN*5.79E-6)
  RCONST(1642) = (1.00D+06)
  RCONST(1643) = (1.02D-11)
  RCONST(1644) = (SUN*5.79E-6+SUN*1.99E-5)
  RCONST(1645) = (1.58D-11)
  RCONST(1646) = (2.90D-12)
  RCONST(1647) = (SUN*3.98E-6)
  RCONST(1648) = (2.91D-13*EXP(1300/TEMP))
  RCONST(1649) = (2.7D-12*EXP(360/TEMP))
! RCONST(1650) = constant rate coefficient
  RCONST(1651) = (2.00D-12*RO2)
  RCONST(1652) = (1.4D-12*EXP(-1860/TEMP)*5.5)
  RCONST(1653) = (7.09D-11)
  RCONST(1654) = (SUN*1.99E-5)
  RCONST(1655) = (3.59D-12)
  RCONST(1656) = (SUN*5.79E-6)
  RCONST(1657) = (1.00D+06)
  RCONST(1658) = (8.76D-13)
  RCONST(1659) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1660) = (2.7D-12*EXP(360/TEMP))
! RCONST(1661) = constant rate coefficient
  RCONST(1662) = (8.00D-13*RO2)
  RCONST(1663) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1664) = (2.7D-12*EXP(360/TEMP))
! RCONST(1665) = constant rate coefficient
  RCONST(1666) = (8.00D-13*RO2)
  RCONST(1667) = (1.53D-12)
  RCONST(1668) = (SUN*5.79E-6)
  RCONST(1669) = (1.05D-10)
  RCONST(1670) = (SUN*5.79E-6)
  RCONST(1671) = (1.00D+06)
  RCONST(1672) = (1.61D-10)
  RCONST(1673) = (SUN*5.79E-6)
  RCONST(1674) = (1.00D+06)
  RCONST(1675) = (9.85D-11)
  RCONST(1676) = (1.30D-10)
  RCONST(1677) = (1.90D-12*EXP(190/TEMP))
  RCONST(1678) = (SUN*5.79E-6)
  RCONST(1679) = (1.00D+06)
  RCONST(1680) = (1.90D-12*EXP(190/TEMP))
  RCONST(1681) = (SUN*5.79E-6)
  RCONST(1682) = (1.00D+06)
  RCONST(1683) = (4.06D-12)
  RCONST(1684) = (SUN*3.98E-6*2)
  RCONST(1685) = (6.66D-11)
  RCONST(1686) = (SUN*2.68E-4)
  RCONST(1687) = (1.4D-12*EXP(-1860/TEMP)*8.5)
  RCONST(1688) = (2.68D-11)
  RCONST(1689) = (SUN*1.99E-5)
  RCONST(1690) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1691) = (2.7D-12*EXP(360/TEMP))
! RCONST(1692) = constant rate coefficient
  RCONST(1693) = (6.00D-13*RO2)
  RCONST(1694) = (1.00D+06)
  RCONST(1695) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1696) = (2.7D-12*EXP(360/TEMP))
! RCONST(1697) = constant rate coefficient
  RCONST(1698) = (9.20D-14*RO2)
  RCONST(1699) = (1.4D-12*EXP(-1860/TEMP)*5.5)
  RCONST(1700) = (2.15D-11)
  RCONST(1701) = (SUN*1.99E-5)
  RCONST(1702) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1703) = (2.7D-12*EXP(360/TEMP))
! RCONST(1704) = constant rate coefficient
  RCONST(1705) = (9.20D-14*RO2*0.70)
  RCONST(1706) = (9.20D-14*RO2*0.30)
  RCONST(1707) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1708) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1709) = (7.5D-12*EXP(290/TEMP))
  RCONST(1710) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1711) = (2.3E-12*1.74)
  RCONST(1712) = (1.00D-11*RO2)
  RCONST(1713) = (3.42D-11)
  RCONST(1714) = (SUN*5.79E-6+SUN*1.99E-5)
  RCONST(1715) = (SUN*5.79E-6)
  RCONST(1716) = (1.00D+06)
  RCONST(1717) = (2.91D-13*EXP(1300/TEMP)*0.625)
  RCONST(1718) = (2.7D-12*EXP(360/TEMP))
! RCONST(1719) = constant rate coefficient
  RCONST(1720) = (8.80D-13*0.2*RO2)
  RCONST(1721) = (8.80D-13*0.6*RO2)
  RCONST(1722) = (8.80D-13*0.2*RO2)
  RCONST(1723) = (6.99D-11)
  RCONST(1724) = (SUN*5.79E-6+SUN*5.50E-5)
  RCONST(1725) = (6.64D-11)
  RCONST(1726) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1727) = (SUN*3.98E-6)
  RCONST(1728) = (SUN*5.79E-6)
  RCONST(1729) = (7.34D-12)
  RCONST(1730) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1731) = (3.74D-12)
  RCONST(1732) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1733) = (2.7D-12*EXP(360/TEMP))
! RCONST(1734) = constant rate coefficient
  RCONST(1735) = (8.00D-13*RO2)
  RCONST(1736) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1737) = (2.7D-12*EXP(360/TEMP))
! RCONST(1738) = constant rate coefficient
  RCONST(1739) = (8.00D-13*RO2)
  RCONST(1740) = (1.53D-12)
  RCONST(1741) = (SUN*5.79E-6)
  RCONST(1742) = (7.54D-11)
  RCONST(1743) = (SUN*5.79E-6)
  RCONST(1744) = (1.00D+06)
  RCONST(1745) = (9.48D-11)
  RCONST(1746) = (SUN*5.79E-6)
  RCONST(1747) = (1.00D+06)
  RCONST(1748) = (9.13D-11)
  RCONST(1749) = (1.90D-12*EXP(190/TEMP))
  RCONST(1750) = (SUN*5.79E-6)
  RCONST(1751) = (1.00D+06)
  RCONST(1752) = (1.90D-12*EXP(190/TEMP))
  RCONST(1753) = (SUN*5.79E-6)
  RCONST(1754) = (1.00D+06)
  RCONST(1755) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1756) = (2.7D-12*EXP(360/TEMP))
! RCONST(1757) = constant rate coefficient
  RCONST(1758) = (9.20D-14*RO2)
  RCONST(1759) = (1.00D+06)
  RCONST(1760) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1761) = (2.7D-12*EXP(360/TEMP))
! RCONST(1762) = constant rate coefficient
  RCONST(1763) = (8.00D-13*RO2)
  RCONST(1764) = (2.91D-13*EXP(1300/TEMP)*0.859)
  RCONST(1765) = (2.7D-12*EXP(360/TEMP))
! RCONST(1766) = constant rate coefficient
  RCONST(1767) = (8.00D-13*RO2)
  RCONST(1768) = (1.53D-12)
  RCONST(1769) = (SUN*5.79E-6)
  RCONST(1770) = (7.54D-11)
  RCONST(1771) = (SUN*5.79E-6)
  RCONST(1772) = (1.00D+06)
  RCONST(1773) = (9.48D-11)
  RCONST(1774) = (SUN*5.79E-6)
  RCONST(1775) = (1.00D+06)
  RCONST(1776) = (9.13D-11)
  RCONST(1777) = (1.90D-12*EXP(190/TEMP))
  RCONST(1778) = (SUN*5.79E-6)
  RCONST(1779) = (1.00D+06)
  RCONST(1780) = (1.90D-12*EXP(190/TEMP))
  RCONST(1781) = (SUN*5.79E-6)
  RCONST(1782) = (1.00D+06)
  RCONST(1783) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1784) = (2.7D-12*EXP(360/TEMP))
! RCONST(1785) = constant rate coefficient
  RCONST(1786) = (8.00D-13*RO2)
  RCONST(1787) = (2.91D-13*EXP(1300/TEMP)*0.770)
  RCONST(1788) = (2.7D-12*EXP(360/TEMP))
! RCONST(1789) = constant rate coefficient
  RCONST(1790) = (8.00D-13*RO2)
  RCONST(1791) = (9.00D-13)
  RCONST(1792) = (SUN*5.79E-6)
  RCONST(1793) = (1.90D-12*EXP(190/TEMP))
  RCONST(1794) = (SUN*5.79E-6)
  RCONST(1795) = (1.00D+06)
  RCONST(1796) = (1.90D-12*EXP(190/TEMP))
  RCONST(1797) = (SUN*5.79E-6)
  RCONST(1798) = (1.00D+06)
  RCONST(1799) = (6.68D-11)
  RCONST(1800) = (SUN*5.79E-6)
  RCONST(1801) = (1.00D+06)
  RCONST(1802) = (1.23D-10)
  RCONST(1803) = (SUN*5.79E-6)
  RCONST(1804) = (1.00D+06)
  RCONST(1805) = (6.07D-11)
  RCONST(1806) = (9.18D-11)
  RCONST(1807) = (8.56D-11)
  RCONST(1808) = (SUN*5.79E-6)
  RCONST(1809) = (1.00D+06)
  RCONST(1810) = (1.42D-10)
  RCONST(1811) = (SUN*5.79E-6)
  RCONST(1812) = (1.00D+06)
  RCONST(1813) = (7.95D-11)
  RCONST(1814) = (1.11D-10)
  RCONST(1815) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1816) = (2.7D-12*EXP(360/TEMP))
! RCONST(1817) = constant rate coefficient
  RCONST(1818) = (8.00D-13*RO2)
  RCONST(1819) = (2.91D-13*EXP(1300/TEMP)*0.820)
  RCONST(1820) = (2.7D-12*EXP(360/TEMP))
! RCONST(1821) = constant rate coefficient
  RCONST(1822) = (8.00D-13*RO2)
  RCONST(1823) = (1.53D-12)
  RCONST(1824) = (SUN*5.79E-6)
  RCONST(1825) = (1.90D-12*EXP(190/TEMP))
  RCONST(1826) = (SUN*5.79E-6)
  RCONST(1827) = (1.00D+06)
  RCONST(1828) = (1.90D-12*EXP(190/TEMP))
  RCONST(1829) = (SUN*5.79E-6)
  RCONST(1830) = (1.00D+06)
  RCONST(1831) = (2.73D-12)
  RCONST(1832) = (6.19D-12)
  RCONST(1833) = (SUN*5.79E-6)
  RCONST(1834) = (1.12D-12)
  RCONST(1835) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1836) = (1.00D+06)
  RCONST(1837) = (GCARR(2.66E-12,0.0E+00,200.0))
  RCONST(1838) = (GCARR(1.14E-12,0.0E+00,200.0))
  RCONST(1839) = (SUN*6.49E-6)
  RCONST(1840) = (GC_HACOHA(2.15E-12,0.0E+00,305.0))
  RCONST(1841) = (GC_HACOHB(2.15E-12,0.0E+00,305.0))
  RCONST(1842) = (SUN*2.88E-6)
  RCONST(1843) = (GCARR(4.00E-13,0.0E+00,0.0))
  RCONST(1844) = (1.88D-11)
  RCONST(1845) = (1.90D-12*EXP(190/TEMP))
  RCONST(1846) = (SUN*3.98E-6)
  RCONST(1847) = (SUN*5.79E-6)
  RCONST(1848) = (1.00D+06)
  RCONST(1849) = (5.86D-12)
  RCONST(1850) = (SUN*3.98E-6)
  RCONST(1851) = (1.72D-11)
  RCONST(1852) = (SUN*5.79E-6)
  RCONST(1853) = (1.00D+06)
  RCONST(1854) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1855) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1856) = (7.5D-12*EXP(290/TEMP))
  RCONST(1857) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1858) = (2.3E-12*1.74)
  RCONST(1859) = (1.00D-11*RO2)
  RCONST(1860) = (1.90D-12*EXP(190/TEMP))
  RCONST(1861) = (SUN*5.79E-6)
  RCONST(1862) = (1.00D+06)
  RCONST(1863) = (1.90D-12*EXP(190/TEMP))
  RCONST(1864) = (SUN*5.79E-6)
  RCONST(1865) = (1.00D+06)
  RCONST(1866) = (1.15D-10)
  RCONST(1867) = (SUN*1.25E-4*2)
  RCONST(1868) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1869) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1870) = (7.5D-12*EXP(290/TEMP))
  RCONST(1871) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1872) = (2.3E-12*1.74)
  RCONST(1873) = (1.00D-11*RO2)
  RCONST(1874) = (1.90D-12*EXP(190/TEMP))
  RCONST(1875) = (2.04D-13)
  RCONST(1876) = (SUN*2.68E-4)
  RCONST(1877) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1878) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1879) = (7.5D-12*EXP(290/TEMP))
  RCONST(1880) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1881) = (2.3E-12*1.74)
  RCONST(1882) = (1.00D-11*RO2)
  RCONST(1883) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(1884) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(1885) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1886) = (7.5D-12*EXP(290/TEMP))
  RCONST(1887) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1888) = (2.3E-12*1.74)
  RCONST(1889) = (1.00D-11*RO2*0.3)
  RCONST(1890) = (1.00D-11*RO2*0.7)
  RCONST(1891) = (1.19D-11)
  RCONST(1892) = (SUN*5.79E-6)
  RCONST(1893) = (SUN*3.98E-6)
  RCONST(1894) = (1.00D+06)
  RCONST(1895) = (3.42D-11)
  RCONST(1896) = (SUN*5.79E-6)
  RCONST(1897) = (SUN*5.79E-6)
  RCONST(1898) = (SUN*1.99E-5)
  RCONST(1899) = (1.00D+06)
  RCONST(1900) = (5.2D-13*EXP(980/TEMP)*0.15)
  RCONST(1901) = (5.2D-13*EXP(980/TEMP)*0.41)
  RCONST(1902) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1903) = (7.5D-12*EXP(290/TEMP))
  RCONST(1904) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1905) = (2.3E-12*1.74)
  RCONST(1906) = (1.00D-11*RO2*0.3)
  RCONST(1907) = (1.00D-11*RO2*0.7)
  RCONST(1908) = (1.69D-11)
  RCONST(1909) = (SUN*5.79E-6)
  RCONST(1910) = (1.00D+06)
  RCONST(1911) = (1.34D-11)
  RCONST(1912) = (1.02D-13)
  RCONST(1913) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1914) = (4.66D-11)
  RCONST(1915) = (SUN*5.79E-6)
  RCONST(1916) = (1.00D+06)
  RCONST(1917) = (2.55D-11)
  RCONST(1918) = (5.68D-12)
  RCONST(1919) = (1.90D-12*EXP(190/TEMP))
  RCONST(1920) = (SUN*5.79E-6)
  RCONST(1921) = (1.00D+06)
  RCONST(1922) = (1.90D-12*EXP(190/TEMP))
  RCONST(1923) = (SUN*5.79E-6)
  RCONST(1924) = (1.00D+06)
  RCONST(1925) = (7.04D-11)
  RCONST(1926) = (SUN*1.25E-4+SUN*2.68E-4)
  RCONST(1927) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1928) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1929) = (7.5D-12*EXP(290/TEMP))
  RCONST(1930) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1931) = (2.3E-12*1.74)
  RCONST(1932) = (1.00D-11*RO2)
  RCONST(1933) = (1.52D-10)
  RCONST(1934) = (SUN*5.79E-6)
  RCONST(1935) = (1.00D+06)
  RCONST(1936) = (1.90D-12*EXP(190/TEMP))
  RCONST(1937) = (SUN*5.79E-6)
  RCONST(1938) = (1.00D+06)
  RCONST(1939) = (1.90D-12*EXP(190/TEMP))
  RCONST(1940) = (SUN*5.79E-6)
  RCONST(1941) = (1.00D+06)
  RCONST(1942) = (8.35D-11)
  RCONST(1943) = (SUN*1.25E-4+SUN*2.68E-4)
  RCONST(1944) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1945) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1946) = (7.5D-12*EXP(290/TEMP))
  RCONST(1947) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1948) = (2.3E-12*1.74)
  RCONST(1949) = (1.0D-11*RO2)
  RCONST(1950) = (1.90D-12*EXP(190/TEMP))
  RCONST(1951) = (SUN*5.79E-6)
  RCONST(1952) = (1.00D+06)
  RCONST(1953) = (1.90D-12*EXP(190/TEMP))
  RCONST(1954) = (SUN*5.79E-6)
  RCONST(1955) = (1.00D+06)
  RCONST(1956) = (1.90D-12*EXP(190/TEMP))
  RCONST(1957) = (7.70D-11)
  RCONST(1958) = (SUN*1.25E-4*2)
  RCONST(1959) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1960) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1961) = (7.5D-12*EXP(290/TEMP))
  RCONST(1962) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1963) = (2.3E-12*1.74)
  RCONST(1964) = (1.00D-11*RO2)
  RCONST(1965) = (9.58D-11)
  RCONST(1966) = (SUN*1.25E-4*2)
  RCONST(1967) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1968) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1969) = (7.5D-12*EXP(290/TEMP))
  RCONST(1970) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1971) = (2.3E-12*1.74)
  RCONST(1972) = (1.00D-11*RO2)
  RCONST(1973) = (1.90D-12*EXP(190/TEMP))
  RCONST(1974) = (SUN*5.79E-6)
  RCONST(1975) = (1.00D+06)
  RCONST(1976) = (1.90D-12*EXP(190/TEMP))
  RCONST(1977) = (SUN*5.79E-6)
  RCONST(1978) = (1.00D+06)
  RCONST(1979) = (3.59D-12)
  RCONST(1980) = (SUN*5.79E-6)
  RCONST(1981) = (1.00D-14)
  RCONST(1982) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1983) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(1984) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(1985) = (7.5D-12*EXP(290/TEMP))
  RCONST(1986) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(1987) = (2.3E-12*1.74)
  RCONST(1988) = (1.00D-11*RO2)
  RCONST(1989) = (9.83D-11)
  RCONST(1990) = (SUN*5.79E-6)
  RCONST(1991) = (9.51D-11)
  RCONST(1992) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(1993) = (1.00D+06)
  RCONST(1994) = (2.91D-13*EXP(1300/TEMP)*0.706)
  RCONST(1995) = (2.7D-12*EXP(360/TEMP))
! RCONST(1996) = constant rate coefficient
  RCONST(1997) = (2.00D-12*RO2)
  RCONST(1998) = (4.16D-12)
  RCONST(1999) = (SUN*5.79E-6+SUN*2.68E-4)
  RCONST(2000) = (5.66D-13)
  RCONST(2001) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2002) = (6.58D-12)
  RCONST(2003) = (9.69D-12)
  RCONST(2004) = (SUN*5.79E-6)
  RCONST(2005) = (6.46D-12)
  RCONST(2006) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2007) = (8.76D-13)
  RCONST(2008) = (4.34D-12)
  RCONST(2009) = (SUN*5.79E-6)
  RCONST(2010) = (7.36D-13)
  RCONST(2011) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2012) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(2013) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(2014) = (7.5D-12*EXP(290/TEMP))
  RCONST(2015) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(2016) = (2.3E-12*1.74)
  RCONST(2017) = (1.00D-11*RO2)
  RCONST(2018) = (2.15D-11)
  RCONST(2019) = (SUN*3.98E-6+SUN*1.25E-4)
  RCONST(2020) = (2.75D-11)
  RCONST(2021) = (SUN*5.79E-6)
  RCONST(2022) = (2.42D-11)
  RCONST(2023) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2024) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(2025) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(2026) = (7.5D-12*EXP(290/TEMP))
  RCONST(2027) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(2028) = (2.3E-12*1.74)
  RCONST(2029) = (1.00D-11*RO2)
  RCONST(2030) = (2.75D-11)
  RCONST(2031) = (SUN*5.79E-6)
  RCONST(2032) = (2.42D-11)
  RCONST(2033) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2034) = (1.00D+06)
  RCONST(2035) = (5.2D-13*EXP(980/TEMP)*0.56)
  RCONST(2036) = (5.2D-13*EXP(980/TEMP)*0.44)
  RCONST(2037) = (7.5D-12*EXP(290/TEMP))
  RCONST(2038) = (PAN(3.28E-28,0.,-6.87,1.125E-11,0.,-1.105,0.3))
  RCONST(2039) = (2.3E-12*1.74)
  RCONST(2040) = (1.00D-11*RO2)
  RCONST(2041) = (8.01D-11)
  RCONST(2042) = (SUN*5.79E-6)
  RCONST(2043) = (7.66D-11)
  RCONST(2044) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2045) = (9.29D-11)
  RCONST(2046) = (SUN*5.79E-6)
  RCONST(2047) = (8.96D-11)
  RCONST(2048) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2049) = (9.82D-11)
  RCONST(2050) = (4.15D-11)
  RCONST(2051) = (SUN*5.79E-6+SUN*1.25E-4)
  RCONST(2052) = (3.82D-11)
  RCONST(2053) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2054) = (9.59D-12)
  RCONST(2055) = (SUN*5.79E-6+SUN*2.68E-4)
  RCONST(2056) = (1.00D+06)
  RCONST(2057) = (1.4D-12*EXP(-1860/TEMP)*4.0)
  RCONST(2058) = (1.27D-11)
  RCONST(2059) = (SUN*3.98E-6)
  RCONST(2060) = (3.60D-11)
  RCONST(2061) = (SUN*5.79E-6+SUN*1.25E-4)
  RCONST(2062) = (3.28D-11)
  RCONST(2063) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2064) = (3.60D-11)
  RCONST(2065) = (SUN*5.79E-6+SUN*1.25E-4)
  RCONST(2066) = (3.28D-11)
  RCONST(2067) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
  RCONST(2068) = (3.06D-11)
  RCONST(2069) = (SUN*5.79E-6+SUN*1.25E-4)
  RCONST(2070) = (2.74D-11)
  RCONST(2071) = (PAN(1.10E-5,-10100.,0.,1.90E+17,-14100.,0.,0.3))
! RCONST(2072) = constant rate coefficient
! RCONST(2073) = constant rate coefficient
! RCONST(2074) = constant rate coefficient
! RCONST(2075) = constant rate coefficient
  RCONST(2076) = (5.0E-4*2.7558E-7)
  RCONST(2077) = (0.0E-9*1.2401E-8)
  RCONST(2078) = (0.0E-9*1.2023E-8)
  RCONST(2079) = (2.0E-5*6.6138E-9)
! RCONST(2080) = constant rate coefficient
! RCONST(2081) = constant rate coefficient
! RCONST(2082) = constant rate coefficient
! RCONST(2083) = constant rate coefficient
! RCONST(2084) = constant rate coefficient
! RCONST(2085) = constant rate coefficient
! RCONST(2086) = constant rate coefficient
! RCONST(2087) = constant rate coefficient
! RCONST(2088) = constant rate coefficient
! RCONST(2089) = constant rate coefficient
! RCONST(2090) = constant rate coefficient
! RCONST(2091) = constant rate coefficient
! RCONST(2092) = constant rate coefficient
! RCONST(2093) = constant rate coefficient
! RCONST(2094) = constant rate coefficient
! RCONST(2095) = constant rate coefficient
! RCONST(2096) = constant rate coefficient
! RCONST(2097) = constant rate coefficient
! RCONST(2098) = constant rate coefficient
! RCONST(2099) = constant rate coefficient
! RCONST(2100) = constant rate coefficient
! RCONST(2101) = constant rate coefficient
! RCONST(2102) = constant rate coefficient
! RCONST(2103) = constant rate coefficient
! RCONST(2104) = constant rate coefficient
! RCONST(2105) = constant rate coefficient
! RCONST(2106) = constant rate coefficient
! RCONST(2107) = constant rate coefficient
! RCONST(2108) = constant rate coefficient
! RCONST(2109) = constant rate coefficient
! RCONST(2110) = constant rate coefficient
! RCONST(2111) = constant rate coefficient
! RCONST(2112) = constant rate coefficient
! RCONST(2113) = constant rate coefficient
! RCONST(2114) = constant rate coefficient
! RCONST(2115) = constant rate coefficient
! RCONST(2116) = constant rate coefficient
! RCONST(2117) = constant rate coefficient
! RCONST(2118) = constant rate coefficient
! RCONST(2119) = constant rate coefficient
! RCONST(2120) = constant rate coefficient
! RCONST(2121) = constant rate coefficient
! RCONST(2122) = constant rate coefficient
! RCONST(2123) = constant rate coefficient
! RCONST(2124) = constant rate coefficient
! RCONST(2125) = constant rate coefficient
! RCONST(2126) = constant rate coefficient
! RCONST(2127) = constant rate coefficient
! RCONST(2128) = constant rate coefficient
! RCONST(2129) = constant rate coefficient
! RCONST(2130) = constant rate coefficient
! RCONST(2131) = constant rate coefficient
! RCONST(2132) = constant rate coefficient
! RCONST(2133) = constant rate coefficient
! RCONST(2134) = constant rate coefficient
! RCONST(2135) = constant rate coefficient
! RCONST(2136) = constant rate coefficient
! RCONST(2137) = constant rate coefficient
! RCONST(2138) = constant rate coefficient
! RCONST(2139) = constant rate coefficient
! RCONST(2140) = constant rate coefficient
! RCONST(2141) = constant rate coefficient
! RCONST(2142) = constant rate coefficient
! RCONST(2143) = constant rate coefficient
! RCONST(2144) = constant rate coefficient
! RCONST(2145) = constant rate coefficient
! RCONST(2146) = constant rate coefficient
! RCONST(2147) = constant rate coefficient
! RCONST(2148) = constant rate coefficient
! RCONST(2149) = constant rate coefficient
! RCONST(2150) = constant rate coefficient
! RCONST(2151) = constant rate coefficient
! RCONST(2152) = constant rate coefficient
! RCONST(2153) = constant rate coefficient
! RCONST(2154) = constant rate coefficient
! RCONST(2155) = constant rate coefficient
! RCONST(2156) = constant rate coefficient
! RCONST(2157) = constant rate coefficient
! RCONST(2158) = constant rate coefficient
! RCONST(2159) = constant rate coefficient
! RCONST(2160) = constant rate coefficient
! RCONST(2161) = constant rate coefficient
! RCONST(2162) = constant rate coefficient
! RCONST(2163) = constant rate coefficient
! RCONST(2164) = constant rate coefficient
! RCONST(2165) = constant rate coefficient
! RCONST(2166) = constant rate coefficient
! RCONST(2167) = constant rate coefficient
! RCONST(2168) = constant rate coefficient
! RCONST(2169) = constant rate coefficient
! RCONST(2170) = constant rate coefficient
! RCONST(2171) = constant rate coefficient
! RCONST(2172) = constant rate coefficient
! RCONST(2173) = constant rate coefficient
! RCONST(2174) = constant rate coefficient
! RCONST(2175) = constant rate coefficient
! RCONST(2176) = constant rate coefficient
! RCONST(2177) = constant rate coefficient
! RCONST(2178) = constant rate coefficient
! RCONST(2179) = constant rate coefficient
! RCONST(2180) = constant rate coefficient
! RCONST(2181) = constant rate coefficient
! RCONST(2182) = constant rate coefficient
! RCONST(2183) = constant rate coefficient
! RCONST(2184) = constant rate coefficient
! RCONST(2185) = constant rate coefficient
! RCONST(2186) = constant rate coefficient
! RCONST(2187) = constant rate coefficient
! RCONST(2188) = constant rate coefficient
! RCONST(2189) = constant rate coefficient
! RCONST(2190) = constant rate coefficient
! RCONST(2191) = constant rate coefficient
! RCONST(2192) = constant rate coefficient
! RCONST(2193) = constant rate coefficient
! RCONST(2194) = constant rate coefficient
! RCONST(2195) = constant rate coefficient
! RCONST(2196) = constant rate coefficient
! RCONST(2197) = constant rate coefficient
! RCONST(2198) = constant rate coefficient
! RCONST(2199) = constant rate coefficient
! RCONST(2200) = constant rate coefficient
! RCONST(2201) = constant rate coefficient
! RCONST(2202) = constant rate coefficient
! RCONST(2203) = constant rate coefficient
! RCONST(2204) = constant rate coefficient
! RCONST(2205) = constant rate coefficient
! RCONST(2206) = constant rate coefficient
! RCONST(2207) = constant rate coefficient
! RCONST(2208) = constant rate coefficient
! RCONST(2209) = constant rate coefficient
! RCONST(2210) = constant rate coefficient
! RCONST(2211) = constant rate coefficient
! RCONST(2212) = constant rate coefficient
! RCONST(2213) = constant rate coefficient
! RCONST(2214) = constant rate coefficient
! RCONST(2215) = constant rate coefficient
! RCONST(2216) = constant rate coefficient
! RCONST(2217) = constant rate coefficient
! RCONST(2218) = constant rate coefficient
! RCONST(2219) = constant rate coefficient
! RCONST(2220) = constant rate coefficient
! RCONST(2221) = constant rate coefficient
! RCONST(2222) = constant rate coefficient
! RCONST(2223) = constant rate coefficient
! RCONST(2224) = constant rate coefficient
! RCONST(2225) = constant rate coefficient
! RCONST(2226) = constant rate coefficient
! RCONST(2227) = constant rate coefficient
! RCONST(2228) = constant rate coefficient
! RCONST(2229) = constant rate coefficient
! RCONST(2230) = constant rate coefficient
! RCONST(2231) = constant rate coefficient
! RCONST(2232) = constant rate coefficient
! RCONST(2233) = constant rate coefficient
! RCONST(2234) = constant rate coefficient
! RCONST(2235) = constant rate coefficient
! RCONST(2236) = constant rate coefficient
! RCONST(2237) = constant rate coefficient
! RCONST(2238) = constant rate coefficient
! RCONST(2239) = constant rate coefficient
! RCONST(2240) = constant rate coefficient
! RCONST(2241) = constant rate coefficient
! RCONST(2242) = constant rate coefficient
! RCONST(2243) = constant rate coefficient
! RCONST(2244) = constant rate coefficient
! RCONST(2245) = constant rate coefficient
! RCONST(2246) = constant rate coefficient
! RCONST(2247) = constant rate coefficient
! RCONST(2248) = constant rate coefficient
! RCONST(2249) = constant rate coefficient
! RCONST(2250) = constant rate coefficient
! RCONST(2251) = constant rate coefficient
! RCONST(2252) = constant rate coefficient
! RCONST(2253) = constant rate coefficient
! RCONST(2254) = constant rate coefficient
! RCONST(2255) = constant rate coefficient
! RCONST(2256) = constant rate coefficient
! RCONST(2257) = constant rate coefficient
! RCONST(2258) = constant rate coefficient
! RCONST(2259) = constant rate coefficient
! RCONST(2260) = constant rate coefficient
! RCONST(2261) = constant rate coefficient
! RCONST(2262) = constant rate coefficient
! RCONST(2263) = constant rate coefficient
! RCONST(2264) = constant rate coefficient
! RCONST(2265) = constant rate coefficient
! RCONST(2266) = constant rate coefficient
! RCONST(2267) = constant rate coefficient
! RCONST(2268) = constant rate coefficient
! RCONST(2269) = constant rate coefficient
! RCONST(2270) = constant rate coefficient
! RCONST(2271) = constant rate coefficient
! RCONST(2272) = constant rate coefficient
! RCONST(2273) = constant rate coefficient
! RCONST(2274) = constant rate coefficient
! RCONST(2275) = constant rate coefficient
! RCONST(2276) = constant rate coefficient
! RCONST(2277) = constant rate coefficient
! RCONST(2278) = constant rate coefficient
! RCONST(2279) = constant rate coefficient
! RCONST(2280) = constant rate coefficient
! RCONST(2281) = constant rate coefficient
! RCONST(2282) = constant rate coefficient
! RCONST(2283) = constant rate coefficient
! RCONST(2284) = constant rate coefficient
! RCONST(2285) = constant rate coefficient
! RCONST(2286) = constant rate coefficient
! RCONST(2287) = constant rate coefficient
! RCONST(2288) = constant rate coefficient
! RCONST(2289) = constant rate coefficient
! RCONST(2290) = constant rate coefficient
! RCONST(2291) = constant rate coefficient
! RCONST(2292) = constant rate coefficient
! RCONST(2293) = constant rate coefficient
! RCONST(2294) = constant rate coefficient
! RCONST(2295) = constant rate coefficient
! RCONST(2296) = constant rate coefficient
! RCONST(2297) = constant rate coefficient
! RCONST(2298) = constant rate coefficient
! RCONST(2299) = constant rate coefficient
! RCONST(2300) = constant rate coefficient
! RCONST(2301) = constant rate coefficient
! RCONST(2302) = constant rate coefficient
! RCONST(2303) = constant rate coefficient
! RCONST(2304) = constant rate coefficient
! RCONST(2305) = constant rate coefficient
! RCONST(2306) = constant rate coefficient
! RCONST(2307) = constant rate coefficient
! RCONST(2308) = constant rate coefficient
! RCONST(2309) = constant rate coefficient
! RCONST(2310) = constant rate coefficient
! RCONST(2311) = constant rate coefficient
! RCONST(2312) = constant rate coefficient
! RCONST(2313) = constant rate coefficient
! RCONST(2314) = constant rate coefficient
! RCONST(2315) = constant rate coefficient
! RCONST(2316) = constant rate coefficient
! RCONST(2317) = constant rate coefficient
! RCONST(2318) = constant rate coefficient
! RCONST(2319) = constant rate coefficient
! RCONST(2320) = constant rate coefficient
! RCONST(2321) = constant rate coefficient
! RCONST(2322) = constant rate coefficient
! RCONST(2323) = constant rate coefficient
! RCONST(2324) = constant rate coefficient
! RCONST(2325) = constant rate coefficient
! RCONST(2326) = constant rate coefficient
! RCONST(2327) = constant rate coefficient
! RCONST(2328) = constant rate coefficient
! RCONST(2329) = constant rate coefficient
! RCONST(2330) = constant rate coefficient
! RCONST(2331) = constant rate coefficient
! RCONST(2332) = constant rate coefficient
! RCONST(2333) = constant rate coefficient
! RCONST(2334) = constant rate coefficient
! RCONST(2335) = constant rate coefficient
! RCONST(2336) = constant rate coefficient
! RCONST(2337) = constant rate coefficient
! RCONST(2338) = constant rate coefficient
! RCONST(2339) = constant rate coefficient
! RCONST(2340) = constant rate coefficient
! RCONST(2341) = constant rate coefficient
! RCONST(2342) = constant rate coefficient
! RCONST(2343) = constant rate coefficient
! RCONST(2344) = constant rate coefficient
! RCONST(2345) = constant rate coefficient
! RCONST(2346) = constant rate coefficient
! RCONST(2347) = constant rate coefficient
! RCONST(2348) = constant rate coefficient
! RCONST(2349) = constant rate coefficient
! RCONST(2350) = constant rate coefficient
! RCONST(2351) = constant rate coefficient
! RCONST(2352) = constant rate coefficient
! RCONST(2353) = constant rate coefficient
! RCONST(2354) = constant rate coefficient
! RCONST(2355) = constant rate coefficient
! RCONST(2356) = constant rate coefficient
! RCONST(2357) = constant rate coefficient
! RCONST(2358) = constant rate coefficient
! RCONST(2359) = constant rate coefficient
! RCONST(2360) = constant rate coefficient
! RCONST(2361) = constant rate coefficient
! RCONST(2362) = constant rate coefficient
! RCONST(2363) = constant rate coefficient
! RCONST(2364) = constant rate coefficient
! RCONST(2365) = constant rate coefficient
! RCONST(2366) = constant rate coefficient
! RCONST(2367) = constant rate coefficient
! RCONST(2368) = constant rate coefficient
! RCONST(2369) = constant rate coefficient
! RCONST(2370) = constant rate coefficient
! RCONST(2371) = constant rate coefficient
! RCONST(2372) = constant rate coefficient
! RCONST(2373) = constant rate coefficient
! RCONST(2374) = constant rate coefficient
! RCONST(2375) = constant rate coefficient
! RCONST(2376) = constant rate coefficient
! RCONST(2377) = constant rate coefficient
! RCONST(2378) = constant rate coefficient
! RCONST(2379) = constant rate coefficient
! RCONST(2380) = constant rate coefficient
! RCONST(2381) = constant rate coefficient
! RCONST(2382) = constant rate coefficient
! RCONST(2383) = constant rate coefficient
! RCONST(2384) = constant rate coefficient
! RCONST(2385) = constant rate coefficient
! RCONST(2386) = constant rate coefficient
! RCONST(2387) = constant rate coefficient
! RCONST(2388) = constant rate coefficient
! RCONST(2389) = constant rate coefficient
! RCONST(2390) = constant rate coefficient
! RCONST(2391) = constant rate coefficient
! RCONST(2392) = constant rate coefficient
! RCONST(2393) = constant rate coefficient
! RCONST(2394) = constant rate coefficient
! RCONST(2395) = constant rate coefficient
! RCONST(2396) = constant rate coefficient
! RCONST(2397) = constant rate coefficient
! RCONST(2398) = constant rate coefficient
! RCONST(2399) = constant rate coefficient
! RCONST(2400) = constant rate coefficient
! RCONST(2401) = constant rate coefficient
! RCONST(2402) = constant rate coefficient
! RCONST(2403) = constant rate coefficient
! RCONST(2404) = constant rate coefficient
! RCONST(2405) = constant rate coefficient
! RCONST(2406) = constant rate coefficient
! RCONST(2407) = constant rate coefficient
! RCONST(2408) = constant rate coefficient
! RCONST(2409) = constant rate coefficient
! RCONST(2410) = constant rate coefficient
! RCONST(2411) = constant rate coefficient
! RCONST(2412) = constant rate coefficient
! RCONST(2413) = constant rate coefficient
! RCONST(2414) = constant rate coefficient
! RCONST(2415) = constant rate coefficient
! RCONST(2416) = constant rate coefficient
! RCONST(2417) = constant rate coefficient
! RCONST(2418) = constant rate coefficient
! RCONST(2419) = constant rate coefficient
! RCONST(2420) = constant rate coefficient
! RCONST(2421) = constant rate coefficient
! RCONST(2422) = constant rate coefficient
! RCONST(2423) = constant rate coefficient
! RCONST(2424) = constant rate coefficient
! RCONST(2425) = constant rate coefficient
! RCONST(2426) = constant rate coefficient
! RCONST(2427) = constant rate coefficient
! RCONST(2428) = constant rate coefficient
! RCONST(2429) = constant rate coefficient
! RCONST(2430) = constant rate coefficient
! RCONST(2431) = constant rate coefficient
! RCONST(2432) = constant rate coefficient
! RCONST(2433) = constant rate coefficient
! RCONST(2434) = constant rate coefficient
! RCONST(2435) = constant rate coefficient
! RCONST(2436) = constant rate coefficient
! RCONST(2437) = constant rate coefficient
! RCONST(2438) = constant rate coefficient
! RCONST(2439) = constant rate coefficient
! RCONST(2440) = constant rate coefficient
! RCONST(2441) = constant rate coefficient
! RCONST(2442) = constant rate coefficient
! RCONST(2443) = constant rate coefficient
! RCONST(2444) = constant rate coefficient
! RCONST(2445) = constant rate coefficient
! RCONST(2446) = constant rate coefficient
! RCONST(2447) = constant rate coefficient
! RCONST(2448) = constant rate coefficient
! RCONST(2449) = constant rate coefficient
! RCONST(2450) = constant rate coefficient
! RCONST(2451) = constant rate coefficient
! RCONST(2452) = constant rate coefficient
! RCONST(2453) = constant rate coefficient
! RCONST(2454) = constant rate coefficient
! RCONST(2455) = constant rate coefficient
! RCONST(2456) = constant rate coefficient
! RCONST(2457) = constant rate coefficient
! RCONST(2458) = constant rate coefficient
! RCONST(2459) = constant rate coefficient
! RCONST(2460) = constant rate coefficient
! RCONST(2461) = constant rate coefficient
! RCONST(2462) = constant rate coefficient
! RCONST(2463) = constant rate coefficient
! RCONST(2464) = constant rate coefficient
! RCONST(2465) = constant rate coefficient
! RCONST(2466) = constant rate coefficient
! RCONST(2467) = constant rate coefficient
! RCONST(2468) = constant rate coefficient
! RCONST(2469) = constant rate coefficient
! RCONST(2470) = constant rate coefficient
! RCONST(2471) = constant rate coefficient
! RCONST(2472) = constant rate coefficient
! RCONST(2473) = constant rate coefficient
! RCONST(2474) = constant rate coefficient
! RCONST(2475) = constant rate coefficient
! RCONST(2476) = constant rate coefficient
! RCONST(2477) = constant rate coefficient
! RCONST(2478) = constant rate coefficient
! RCONST(2479) = constant rate coefficient
! RCONST(2480) = constant rate coefficient
! RCONST(2481) = constant rate coefficient
! RCONST(2482) = constant rate coefficient
! RCONST(2483) = constant rate coefficient
! RCONST(2484) = constant rate coefficient
! RCONST(2485) = constant rate coefficient
! RCONST(2486) = constant rate coefficient
! RCONST(2487) = constant rate coefficient
! RCONST(2488) = constant rate coefficient
! RCONST(2489) = constant rate coefficient
! RCONST(2490) = constant rate coefficient
! RCONST(2491) = constant rate coefficient
! RCONST(2492) = constant rate coefficient
! RCONST(2493) = constant rate coefficient
! RCONST(2494) = constant rate coefficient
! RCONST(2495) = constant rate coefficient
! RCONST(2496) = constant rate coefficient
! RCONST(2497) = constant rate coefficient
! RCONST(2498) = constant rate coefficient
! RCONST(2499) = constant rate coefficient
! RCONST(2500) = constant rate coefficient
! RCONST(2501) = constant rate coefficient
! RCONST(2502) = constant rate coefficient
! RCONST(2503) = constant rate coefficient
! RCONST(2504) = constant rate coefficient
! RCONST(2505) = constant rate coefficient
! RCONST(2506) = constant rate coefficient
! RCONST(2507) = constant rate coefficient
! RCONST(2508) = constant rate coefficient
! RCONST(2509) = constant rate coefficient
! RCONST(2510) = constant rate coefficient
! RCONST(2511) = constant rate coefficient
! RCONST(2512) = constant rate coefficient
! RCONST(2513) = constant rate coefficient
! RCONST(2514) = constant rate coefficient
! RCONST(2515) = constant rate coefficient
! RCONST(2516) = constant rate coefficient
! RCONST(2517) = constant rate coefficient
! RCONST(2518) = constant rate coefficient
! RCONST(2519) = constant rate coefficient
! RCONST(2520) = constant rate coefficient
! RCONST(2521) = constant rate coefficient
! RCONST(2522) = constant rate coefficient
! RCONST(2523) = constant rate coefficient
! RCONST(2524) = constant rate coefficient
! RCONST(2525) = constant rate coefficient
! RCONST(2526) = constant rate coefficient
! RCONST(2527) = constant rate coefficient
! RCONST(2528) = constant rate coefficient
! RCONST(2529) = constant rate coefficient
! RCONST(2530) = constant rate coefficient
! RCONST(2531) = constant rate coefficient
! RCONST(2532) = constant rate coefficient
! RCONST(2533) = constant rate coefficient
! RCONST(2534) = constant rate coefficient
! RCONST(2535) = constant rate coefficient
! RCONST(2536) = constant rate coefficient
! RCONST(2537) = constant rate coefficient
! RCONST(2538) = constant rate coefficient
! RCONST(2539) = constant rate coefficient
! RCONST(2540) = constant rate coefficient
! RCONST(2541) = constant rate coefficient
! RCONST(2542) = constant rate coefficient
! RCONST(2543) = constant rate coefficient
! RCONST(2544) = constant rate coefficient
! RCONST(2545) = constant rate coefficient
! RCONST(2546) = constant rate coefficient
! RCONST(2547) = constant rate coefficient
! RCONST(2548) = constant rate coefficient
! RCONST(2549) = constant rate coefficient
! RCONST(2550) = constant rate coefficient
! RCONST(2551) = constant rate coefficient
! RCONST(2552) = constant rate coefficient
! RCONST(2553) = constant rate coefficient
! RCONST(2554) = constant rate coefficient
! RCONST(2555) = constant rate coefficient
! RCONST(2556) = constant rate coefficient
! RCONST(2557) = constant rate coefficient
! RCONST(2558) = constant rate coefficient
! RCONST(2559) = constant rate coefficient
! RCONST(2560) = constant rate coefficient
! RCONST(2561) = constant rate coefficient
! RCONST(2562) = constant rate coefficient
! RCONST(2563) = constant rate coefficient
! RCONST(2564) = constant rate coefficient
! RCONST(2565) = constant rate coefficient
! RCONST(2566) = constant rate coefficient
! RCONST(2567) = constant rate coefficient
! RCONST(2568) = constant rate coefficient
! RCONST(2569) = constant rate coefficient
! RCONST(2570) = constant rate coefficient
! RCONST(2571) = constant rate coefficient
! RCONST(2572) = constant rate coefficient
! RCONST(2573) = constant rate coefficient
! RCONST(2574) = constant rate coefficient
! RCONST(2575) = constant rate coefficient
! RCONST(2576) = constant rate coefficient
! RCONST(2577) = constant rate coefficient
! RCONST(2578) = constant rate coefficient
! RCONST(2579) = constant rate coefficient
! RCONST(2580) = constant rate coefficient
! RCONST(2581) = constant rate coefficient
! RCONST(2582) = constant rate coefficient
! RCONST(2583) = constant rate coefficient
! RCONST(2584) = constant rate coefficient
! RCONST(2585) = constant rate coefficient
! RCONST(2586) = constant rate coefficient
! RCONST(2587) = constant rate coefficient
! RCONST(2588) = constant rate coefficient
! RCONST(2589) = constant rate coefficient
! RCONST(2590) = constant rate coefficient
! RCONST(2591) = constant rate coefficient
! RCONST(2592) = constant rate coefficient
! RCONST(2593) = constant rate coefficient
! RCONST(2594) = constant rate coefficient
! RCONST(2595) = constant rate coefficient
! RCONST(2596) = constant rate coefficient
! RCONST(2597) = constant rate coefficient
! RCONST(2598) = constant rate coefficient
! RCONST(2599) = constant rate coefficient
! RCONST(2600) = constant rate coefficient
! RCONST(2601) = constant rate coefficient
! RCONST(2602) = constant rate coefficient
! RCONST(2603) = constant rate coefficient
! RCONST(2604) = constant rate coefficient
! RCONST(2605) = constant rate coefficient
! RCONST(2606) = constant rate coefficient
! RCONST(2607) = constant rate coefficient
! RCONST(2608) = constant rate coefficient
! RCONST(2609) = constant rate coefficient
! RCONST(2610) = constant rate coefficient
! RCONST(2611) = constant rate coefficient
! RCONST(2612) = constant rate coefficient
! RCONST(2613) = constant rate coefficient
! RCONST(2614) = constant rate coefficient
! RCONST(2615) = constant rate coefficient
! RCONST(2616) = constant rate coefficient
! RCONST(2617) = constant rate coefficient
! RCONST(2618) = constant rate coefficient
! RCONST(2619) = constant rate coefficient
! RCONST(2620) = constant rate coefficient
! RCONST(2621) = constant rate coefficient
! RCONST(2622) = constant rate coefficient
! RCONST(2623) = constant rate coefficient
! RCONST(2624) = constant rate coefficient
! RCONST(2625) = constant rate coefficient
! RCONST(2626) = constant rate coefficient
! RCONST(2627) = constant rate coefficient
! RCONST(2628) = constant rate coefficient
! RCONST(2629) = constant rate coefficient
! RCONST(2630) = constant rate coefficient
! RCONST(2631) = constant rate coefficient
! RCONST(2632) = constant rate coefficient
! RCONST(2633) = constant rate coefficient
! RCONST(2634) = constant rate coefficient
! RCONST(2635) = constant rate coefficient
! RCONST(2636) = constant rate coefficient
! RCONST(2637) = constant rate coefficient
! RCONST(2638) = constant rate coefficient
! RCONST(2639) = constant rate coefficient
! RCONST(2640) = constant rate coefficient
! RCONST(2641) = constant rate coefficient
! RCONST(2642) = constant rate coefficient
! RCONST(2643) = constant rate coefficient
! RCONST(2644) = constant rate coefficient
! RCONST(2645) = constant rate coefficient
! RCONST(2646) = constant rate coefficient
! RCONST(2647) = constant rate coefficient
! RCONST(2648) = constant rate coefficient
! RCONST(2649) = constant rate coefficient
! RCONST(2650) = constant rate coefficient
! RCONST(2651) = constant rate coefficient
! RCONST(2652) = constant rate coefficient
! RCONST(2653) = constant rate coefficient
! RCONST(2654) = constant rate coefficient
! RCONST(2655) = constant rate coefficient
! RCONST(2656) = constant rate coefficient
! RCONST(2657) = constant rate coefficient
! RCONST(2658) = constant rate coefficient
! RCONST(2659) = constant rate coefficient
! RCONST(2660) = constant rate coefficient
! RCONST(2661) = constant rate coefficient
! RCONST(2662) = constant rate coefficient
! RCONST(2663) = constant rate coefficient
! RCONST(2664) = constant rate coefficient
! RCONST(2665) = constant rate coefficient
! RCONST(2666) = constant rate coefficient
! RCONST(2667) = constant rate coefficient
! RCONST(2668) = constant rate coefficient
! RCONST(2669) = constant rate coefficient
! RCONST(2670) = constant rate coefficient
! RCONST(2671) = constant rate coefficient
! RCONST(2672) = constant rate coefficient
! RCONST(2673) = constant rate coefficient
! RCONST(2674) = constant rate coefficient
! RCONST(2675) = constant rate coefficient
! RCONST(2676) = constant rate coefficient
! RCONST(2677) = constant rate coefficient
! RCONST(2678) = constant rate coefficient
! RCONST(2679) = constant rate coefficient
! RCONST(2680) = constant rate coefficient
! RCONST(2681) = constant rate coefficient
! RCONST(2682) = constant rate coefficient
! RCONST(2683) = constant rate coefficient
! RCONST(2684) = constant rate coefficient
! RCONST(2685) = constant rate coefficient
! RCONST(2686) = constant rate coefficient
! RCONST(2687) = constant rate coefficient
! RCONST(2688) = constant rate coefficient
! RCONST(2689) = constant rate coefficient
! RCONST(2690) = constant rate coefficient
! RCONST(2691) = constant rate coefficient
! RCONST(2692) = constant rate coefficient
! RCONST(2693) = constant rate coefficient
! RCONST(2694) = constant rate coefficient
! RCONST(2695) = constant rate coefficient
! RCONST(2696) = constant rate coefficient
! RCONST(2697) = constant rate coefficient
! RCONST(2698) = constant rate coefficient
! RCONST(2699) = constant rate coefficient
! RCONST(2700) = constant rate coefficient
! RCONST(2701) = constant rate coefficient
! RCONST(2702) = constant rate coefficient
! RCONST(2703) = constant rate coefficient
! RCONST(2704) = constant rate coefficient
! RCONST(2705) = constant rate coefficient
! RCONST(2706) = constant rate coefficient
! RCONST(2707) = constant rate coefficient
! RCONST(2708) = constant rate coefficient
! RCONST(2709) = constant rate coefficient
! RCONST(2710) = constant rate coefficient
! RCONST(2711) = constant rate coefficient
! RCONST(2712) = constant rate coefficient
! RCONST(2713) = constant rate coefficient
! RCONST(2714) = constant rate coefficient
! RCONST(2715) = constant rate coefficient
! RCONST(2716) = constant rate coefficient
! RCONST(2717) = constant rate coefficient
! RCONST(2718) = constant rate coefficient
! RCONST(2719) = constant rate coefficient
! RCONST(2720) = constant rate coefficient
! RCONST(2721) = constant rate coefficient
! RCONST(2722) = constant rate coefficient
! RCONST(2723) = constant rate coefficient
! RCONST(2724) = constant rate coefficient
! RCONST(2725) = constant rate coefficient
! RCONST(2726) = constant rate coefficient
! RCONST(2727) = constant rate coefficient
! RCONST(2728) = constant rate coefficient
! RCONST(2729) = constant rate coefficient
! RCONST(2730) = constant rate coefficient
! RCONST(2731) = constant rate coefficient
! RCONST(2732) = constant rate coefficient
! RCONST(2733) = constant rate coefficient
! RCONST(2734) = constant rate coefficient
! RCONST(2735) = constant rate coefficient
! RCONST(2736) = constant rate coefficient
! RCONST(2737) = constant rate coefficient
! RCONST(2738) = constant rate coefficient
! RCONST(2739) = constant rate coefficient
! RCONST(2740) = constant rate coefficient
! RCONST(2741) = constant rate coefficient
! RCONST(2742) = constant rate coefficient
! RCONST(2743) = constant rate coefficient
! RCONST(2744) = constant rate coefficient
! RCONST(2745) = constant rate coefficient
! RCONST(2746) = constant rate coefficient
! RCONST(2747) = constant rate coefficient
! RCONST(2748) = constant rate coefficient
! RCONST(2749) = constant rate coefficient
! RCONST(2750) = constant rate coefficient
! RCONST(2751) = constant rate coefficient
! RCONST(2752) = constant rate coefficient
! RCONST(2753) = constant rate coefficient
! RCONST(2754) = constant rate coefficient
! RCONST(2755) = constant rate coefficient
! RCONST(2756) = constant rate coefficient
! RCONST(2757) = constant rate coefficient
! RCONST(2758) = constant rate coefficient
! RCONST(2759) = constant rate coefficient
! RCONST(2760) = constant rate coefficient
! RCONST(2761) = constant rate coefficient
! RCONST(2762) = constant rate coefficient
! RCONST(2763) = constant rate coefficient
! RCONST(2764) = constant rate coefficient
! RCONST(2765) = constant rate coefficient
! RCONST(2766) = constant rate coefficient
! RCONST(2767) = constant rate coefficient
! RCONST(2768) = constant rate coefficient
! RCONST(2769) = constant rate coefficient
! RCONST(2770) = constant rate coefficient
! RCONST(2771) = constant rate coefficient
! RCONST(2772) = constant rate coefficient
! RCONST(2773) = constant rate coefficient
! RCONST(2774) = constant rate coefficient
! RCONST(2775) = constant rate coefficient
! RCONST(2776) = constant rate coefficient
! RCONST(2777) = constant rate coefficient
! RCONST(2778) = constant rate coefficient
! RCONST(2779) = constant rate coefficient
! RCONST(2780) = constant rate coefficient
! RCONST(2781) = constant rate coefficient
! RCONST(2782) = constant rate coefficient
! RCONST(2783) = constant rate coefficient
! RCONST(2784) = constant rate coefficient
! RCONST(2785) = constant rate coefficient
! RCONST(2786) = constant rate coefficient
! RCONST(2787) = constant rate coefficient
! RCONST(2788) = constant rate coefficient
! RCONST(2789) = constant rate coefficient
! RCONST(2790) = constant rate coefficient
! RCONST(2791) = constant rate coefficient
! RCONST(2792) = constant rate coefficient
! RCONST(2793) = constant rate coefficient
! RCONST(2794) = constant rate coefficient
! RCONST(2795) = constant rate coefficient
! RCONST(2796) = constant rate coefficient
! RCONST(2797) = constant rate coefficient
! RCONST(2798) = constant rate coefficient
! RCONST(2799) = constant rate coefficient
! RCONST(2800) = constant rate coefficient
! RCONST(2801) = constant rate coefficient
! RCONST(2802) = constant rate coefficient
! RCONST(2803) = constant rate coefficient
! RCONST(2804) = constant rate coefficient
! RCONST(2805) = constant rate coefficient
! RCONST(2806) = constant rate coefficient
! RCONST(2807) = constant rate coefficient
! RCONST(2808) = constant rate coefficient
! RCONST(2809) = constant rate coefficient
! RCONST(2810) = constant rate coefficient
! RCONST(2811) = constant rate coefficient
! RCONST(2812) = constant rate coefficient
      
END SUBROUTINE Update_RCONST

! End of Update_RCONST function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~


! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! Update_PHOTO - function to update photolytical rate constants
!   Arguments :
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

SUBROUTINE Update_PHOTO ( )


   USE aromatics_kpp_Global

      
END SUBROUTINE Update_PHOTO

! End of Update_PHOTO function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



END MODULE aromatics_kpp_Rates

